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摘要 


本 文 首先 简要 地 回顾 了 极 小 曲面 问题 ( Plateau. 问题 ) 的 产生 和 沿革 ， 综 述 了 目前 在 
CAGD 领域 内 Bézier 极 小 曲面 造型 的 主要 方法 ， 如 Dirichlet 泛 函 方法 ， 离 散 的 Mask 方法 
等 . 这 些 方法 是 基于 以 下 结论 : 等 温 正则 参数 曲面 是 极 小 曲面 的 充 要 条 件 是 其 为 调和 曲面 
虽然 这 些 研究 取得 了 一 些 很 好 的 结果 ， 但 是 由 于 Bézier 曲面 不 能 精确 表示 除 抛物 面 外 的 二 
次 曲面 , 只 能 给 出 近似 表示 . 而 有 理 Bézier 曲面 能 克服 这 一 缺陷 ， 所 以 研究 有 理 Bézier Hill 
面 造型 具有 实际 意义 . 基于 此 本 文 进一步 指出 了 研究 空间 闭 曲 线 所 围 成 的 
Plateau -Rational-Bézier 曲面 造型 的 必要 性 和 合理 性 . 本 文 基于 分 而 治之 的 思想 ， 利 用 有 
限 差分 法 来 研究 有 理 Bézier 极 小 曲面 造型 问题 .本 文 的 主要 贡献 与 创新 点 可 以 概括 为 以 下 
JL: 

1. TE CAGD 领域 对 于 有 理 Bézier 极 小 曲面 造型 问题 进行 了 成 功 的 尝试 ,使 得 统一 研 

究 任 意 次 的 有 理 Bézier 极 小 曲面 造型 成 为 可 能 ,这 将 促进 CAD 行业 能 够 在 NURBS 
系统 中 计算 与 绘制 极 小 曲面 ,有 望 对 建筑 、 机 械 等 工程 实际 产生 深远 影响 . 

2. 将 Bézier 调和 曲面 的 研究 推广 到 有 理 形式 . 由 于 有 理 Bézier 曲面 的 复杂 形式 ， 本 
文中 没有 直接 应 用 Bézier 调和 曲面 的 求解 模式 来 做 ， 而 是 采用 了 分 而 治之 的 思想 . 
利用 有 限 差分 法 将 有 理 Bézier 调和 曲面 的 造型 问题 转化 为 求解 一 个 线性 方程 组 ， 
其 可 应 用 到 任意 阶 的 有 理 Bézier 调和 曲面 造型 问题 ， 其 求解 过 程 与 低 阶 情形 一 样 
都 是 求解 一 个 线性 方程 组 

3. 把 上 述 给 定 控制 边界 求解 有 理 Bézier 调和 曲面 的 想法 加 以 编程 实现 ， 通 过 实例 加 
以 验证 ， 说 明 用 有 限 差 分 法 求解 得 到 的 一 类 有 理 Bezier 调和 曲面 最 接近 于 有 理 
Bézier 极 小 曲面 . 另 一 面 ， 为 了 从 客观 上 对 曲面 加 以 细致 分 析 比 较 ， 本 文 给 出 了 其 
曲面 的 绝对 平均 曲率 图 . 


E 
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Abstract 


To begin with, we briefly retrospect the birth and the evolution of the minimal surface 
problem (Plateau Problem) with prescribed border, and in addition we comprehensively introduce 
some important methods, for example, Dirichlet functional method and discrete Mask method and 
so on, of recent researches on Bézier minimal surface modeling in the field of CAGD. These 
methods are based on the following conclusion that for isothermal regular parametric surface, it is 
a minimal surface if and only if this is a harmonic surface. Although these studies have some good 
results, Bézier surface can not be accurately but approximatively expressed for quadric in addition 
to paraboloid. It is well-known that rational Bézier surface can be accurately expressed for not 
only paraboloid but also for quadric. Furthermore, we point out that it is necessary and reasonable 
for us to study Plateau -Rational-Bézier surface modeling surrounded by the closed curve 
problem . Based on the idea of “dividing and conquering”, we apply the method of finite 
difference method to study the problem of rational Bézier minimal surface modeling. The most 
important contributions and innovations of this thesis can be summed up as follows: 

1. We successfully attempt to solve the rational Bézier surfaces modeling problem, which make 
it possible for us to study arbitrary harmonic rational Bézier surfaces modeling. Hence, it will 
promote the field of CAD to use the NURBS system to draw and compute minimal surface, 
which will have a profound influence on some fields of engineering, such as construction and 
mechanism. 

2. We extend the geometric design on the surface of Bézier harmonic to the rational form. 
Owing to the complicated form of rational Bézier surfaces, we don’t directly apply the 
solution to harmonic Bézier surface modeling, but we can adopt the dividing and conquering 
approach. We transform the problem of rational Bézier surfaces modeling into solving linear 
equations. This method can be applied to arbitrary order of the rational Bézier surfaces 
modeling problem, which is the same with low-level case by solving a linear equations. 

3. We implement the algorithm of solving the harmonic rational Bézier surface for given control 
borders and verify this algorithm by several examples in order to explain that a class of 
rational Bézier harmonic surface ,obtained by finite difference method, is the closest to the 
rational Bézier minimal surface. On the other hand, in order to objectively analysis and 
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compare the surfaces, we give their absolute mean curvature maps. 


Keywords: minimal surface, Rational Bézier surface, Plateau-Rational-Bézier problem, 


Dirichlet energy function, Harmonic Surfaces, Isothermal parameter surface, 


Finite difference method. 
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第 一 章 绪论 


$ 1.1 极 小 曲面 的 提出 及 研究 


至 一 个 藏 于 罗浮 襄 的 古 
液 ， 并 带 有 一 个 连带 手相 


的 简单 的 金 


一 点 ， 类 似 圆 盘 形 状 的 肥 拒 膜 就 形成 了 ， 
的 几何 形状 的 架子 以 替代 
涉 效应 而 导致 彩虹 色 ， 但 


杂 


的 金属 丝 做 成 更 
的 波动 性 ， 肥 皂 膜 会 因 干 


及 顿 花瓶 上 也 会 有 


在 我 们 的 孩童 时 代 ， 很 平常 的 游戏 就 是 吹 肥 是 泡 . 吹 肥 气泡 是 一 种 古老 的 娱乐 方式 ， 甚 
吹 肥 气 泡 的 儿童 画像 . 今 可 购 得 
遇 丝 架 或 塑料 架子 ， 把 这 种 架子 浸入 浴 液 ， 然 后 沾 起 


特 


制 的 泡沫 浴 


向 架子 吹 气 ， 就 产 


AAT, N 


生 肥 气泡 . WE 
I 会 出 现 更 加 神奇 的 效果 .由 


32h 


瞬息 后 


了 .肥皂 膜 实验 中 出 现 的 诱 人 注 膜 , 显然 是 处 于 稳定 六 
为 极 小 势能 
面 的 面积 比 跨 在 同一 架子 上 而 形成 的 任何 


IRE 
面积 最 小 的 曲面 . 极 小 


[Langrange,1788]， 称 殿 


F 衡 状态 ， 


! 于 注 膜 大 注 ， 七 彩色 便 看 不 到 
那么 按照 约翰 . 伯 努 利 的 虚 功 


Labs 


势能 正比 于 面积 ， 


AL re Hk Be H 


些 好 的 性 质 ? 这 种 


数学 模型 的 奇 漠 形状 直观 可 见 ， 也 使 得 逢 
[Stefan Hildebrandt,2004]. 那么 我 们 从 数学 上 应 当 怎 样 来 
面 在 工程 造型 中 有 应 用 价值 吗 ? 数学 家 们 对 


数学 家 们 称 上 


多 数学 家 对 


述 这 样 


造型 的 数学 曲面 是 
邻近 曲面 的 面积 都 小 . 
极 小 原理 十 分 感 兴趣 
的 曲面 呢 ? 这 种 曲面 有 哪 
此 的 研究 状况 又 如 何 


由 肥 电 


? 


述 薄 膜 曲 面 为 极 小 曲面 ， 它 是 所 有 给 定 特定 边 


界 中 面积 最 小 的 曲 


面 . 极 小 曲面 简单 说 就 是 平均 曲率 为 零 的 曲面 .在 这 


平均 曲率 (mean curvature) 


是 微分 儿 何 中 一 个 “外 在 的 ” 玖 曲 测量 标准 , 局 部 地 描述 了 一 个 曲面 做 入 周围 空间 (〈 比 
如 二 维 曲面 嵌 入 三 维 欧 几 里 得 空间 ) 的 曲率 .这 个 概念 由 索 菲 . 热 尔 曼 在 她 的 著作 《 弹 
性 理论 》 中 最 先 引 入 [Germain ,1831]. 平均 由 率 定 义 为 吾 = (eie) ,其 HE eue, K 
示 两 个 主 曲率 [C. Spivak,1999]. 

给 定 空间 中 的 一 条 可 求 长 (可 测 ) Jordan 闭 曲 线 (简单 闭 曲 线 ,或 者 称 无 重点 的 连 


忆 曲 线 ), 寻 找 以 该 闭 曲 线 为 边界 的 


其 面积 


可 展 曲 面 . 除 平 面 外 ， 


是 正 螺 面 [do Carmo,1976]. 


达到 最 小 一 一 这 一 问题 称 为 极 小 曲面 问 


面 ,使 得 该 


旋转 极 小 曲面 都 是 悬 链 面 [do Carmo,1976]， 


曲面 在 所 有 具有 相同 边界 的 曲面 中 ， 
题 [Nitsche,1988]. 平 面 是 仅 有 的 极 小 


直 纹 极 小 曲面 都 
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极 小 曲面 在 物理 上 解释 为 “肥皂 泡 现象 ”， 然 而 当 把 圆 环 扭曲 但 仍 闭 合 时 , 极 小 曲面 的 


形状 就 不 是 很 明显 了 [do Carmo,1976] .这 是 从 物理 层面 上 对 极 小 曲面 的 揭示 , 很 多 问题 都 


是 


于 奇异 的 现象 引发 的 ， 很 多 问题 都 是 从 直观 到 抽象 的 认识 . 
极 小 曲面 在 工程 领域 有 着 广泛 应 用 ， 在 符合 给 定 边界 条 件 下 求 出 相应 极 小 曲面 是 有 很 
大 意义 的 ， 既 满足 条 件 又 可 以 节约 材料 ， 在 如 今 资源 短缺 的 时 代 是 十 分 关键 的 . 极 小 曲面 
问题 是 微分 几何 领域 一 个 古老 而 活跃 的 主题 [Nitsche,1988]. 对 这 个 问题 研究 最 早 可 以 奶 漳 


到 17 世纪 .， 1744 年 ， 瑞 士 数 学 家 Euler 导出 了 变 分 法 的 Euler 方程 , 发 现 了 某 些 极 小 曲面 


[Euler,1744]. Euler 的 方法 后 来 又 为 法 国 数学 家 Lagrange 所 发 展 ，Lagrange 首先 将 变 分 


法 置 于 分 析 的 基础 上 ， 他 充分 应 用 变 分 法 来 构筑 其 分 析 力 学 体系 ， 并 系统 地 研究 了 变 分 法 
及 其 在 力学 上 的 应 用 ， 全 部 力学 被 他 化 归 为 一 个 统一 的 变 分 原理 一 一 虚 功 原理 


[Lagrange ,1788]. 值得 指出 的 是 ，Lagrange 在 1760 年 就 已 经 用 变 分 原理 推导 出 极 小 曲 


面 应 该 满足 的 偏 微分 方程 ， 称 为 Euler — Lagrange 方程 [GIL BARG,1983]， 该 方程 是 一 个 


拟 线性 二 阶 椭圆 偏 微分 方程 . 极 小 曲面 的 几何 意义 为 曲面 的 平均 曲率 等 于 0, 所 以 极 小 曲面 
就 是 平均 曲率 处 处 为 0 的 曲面 ， 除 了 在 有 限 个 点 以 外 ， 极 小 曲面 的 Grass 曲率 恒 为 负 ， 所 
以 极 小 曲面 具有 许多 优良 的 性 质 ， 其 结构 稳定 ， 自 然 光 顺 ， 在 建筑 设计 、 飞 机 轮船 制造 、 
分 子 化 学 、 晶 体 学 等 各 个 领域 有 很 重要 的 应 用 , 例如 ， 用 极 小 曲面 作 房 项 曲面 ， 结 构 很 稳 
定 ， 美 观 大 方 ， 特 别 有 不 积 水 的 特点 [Man,Wang,2001]. 

早期 就 有 很 多 研究 极 小 曲面 的 数学 家 都 致力 于 寻找 该 偏 微分 方程 闭 形式 的 解 ， 当 然 也 


a 
b 


取得 了 很 多 成 果 : 如 Monge 等 人 在 18 世纪 末期 建立 了 所 谓 的 Monge-Legendre 表示 


[Monge,1794], 19 世纪 早期 Scherk 给 出 了 一 类 极 小 曲面 即 Scherk 曲面 ， 悬 链 面 和 正 螺 面 


都 是 它 的 特例 [do Carmo,1976] ， 这 也 是 当时 知道 的 仅 有 的 极 小 曲面 .到 19 世纪 中 后 期 有 


很 多 杰出 的 数学 家 包括 Catalan, Weierstrass, Rieman,Enneper,Schwarz, Lie, 


Ribacour,Beltrami 等 都 投入 到 极 小 曲面 的 研究 中 来 ， 得 到 了 极 小 曲面 的 


Werestrass-Enneper 表示 [do Carmo,1976] . Enneper 给 出 了 第 一 个 多 项 式 参数 极 小 曲 


面 即 著名 的 Enneper 曲面 [do Carmo,1976]. 另外 ，S.Lie 和 A.Ribacour 首先 使 用 复 变 


数 几 何 的 方法 研究 极 小 曲面 ， 实 曲面 和 虚 曲 面 的 结合 对 极 小 曲面 理论 的 发 展 具 有 很 深远 的 
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意义 ， 使 得 很 多 问题 可 以 统 


分 方程 解 的 存在 性 问题 .此 时 ， 他 们 的 解 不 能 排除 曲面 可 能 存在 奇 点 ( 即 分 文 点 ) 接着 


现 分 支点 [Osserman,1986]. 然而 即使 是 现在 也 不 能 


Osserman 进一步 证 明了 极 4 


处 理 . 他 们 首次 比较 成 功 地 从 理论 上 解决 了 极 小 曲面 的 偏 微 


Hit 18] AS H 


说 这 一 问题 得 到 了 很 彻底 的 解决 ， 极 小 曲面 的 任何 基于 几何 还 是 物理 上 的 研究 成 


何 学 有 着 很 重要 的 页 献 . 


的 问题 . 如 果 对 极 小 曲面 做 理 


很 多 数学 家 对 极 小 曲面 这 一 特殊 的 曲面 十 分 感 兴 


果 都 对 几 


1 此 提出 了 许多 有 趣 而 


富有 挑战 


想 化 的 物理 实验 [do Carmo,1976] , 这 个 实验 可 以 描述 如 下 : 


把 一 个 金属 丝 环 混入 一 倪 肥 虹 水 中 ， 再 把 它 拿 出 来 后 上 面 便 张 着 一 层 内 着 彩虹 的 平 圆 盘 形 


式 的 薄膜 . 任何 别 的 以 此 圆 环 为 界 的 曲面 ， 不 论 是 露出 皱 袜 的 ， 还 是 强烈 凸 出 的 


更 大 的 面积 . 平 圆 盘 是 张 在 贺 


周 上 的 旺 液 膜 能 表现 出 的 最 小 曲面 ， 它 对 应 于 由 圆 


， 都 会 有 


周 决定 的 


极 小 曲面 , 比利时 盲人 物理 学 家 J.A.F . 普 拉 带 奥 通 过 不 断 地 实验 和 研究 发 现 了 很 多 重要 


在 数学 上 我 们 称 它 为 极 小 曲面 ，Plateau 所 做 的 时 
杂 的 .如 将 圆 环 浸入 同样 的 肥皂 液 两 次 ， 就 会 显示 出 两 种 不 同 的 形式 . 另外 ， 当 将 


力学 》[Plateau,1873]: 通过 实验 断定 ， 神 奇 的 自然 


的 圆 环 上 , 在 忽略 皂 液 膜 重力 的 情况 下 ， 


果 ,， 他 花 了 多 年 的 时 间 做 这 样 的 液体 薄膜 实验 ， 思 索 归 纳 其 明显 的 形状 ， 并 在 1873 年 出 
版 了 著作 《实验 和 理论 流体 前 
到 一 种 方式 来 形成 皂 液 膜 , 使 其 张 在 不 论 如 何 扭 
于 表面 张力 的 作用 ,在 固定 边界 曲线 上 自然 张 成 的 曲面 是 所 有 可 


界 总 能 找 


能 的 曲面 中 面积 最 小 的 ， 


液 膜 实验 实际 上 从 物理 角度 解释 是 很 复 


定 角 度 后 浸入 溶液 中 多 次 ， 每 次 都 会 产生 不 同 的 蛙 液 膜 形 式 . 哲学 家 赫 拉克 利 特 说 :“ 一 


个 人 不 可 能 两 次 跨 过 同一 条 河流 ”. 这 虽然 只 是 一 个 哲学 家 的 观点 ， 但 是 的 确 是 


多 简单 的 物理 现象 ， 在 理论 上 是 很 难 揭示 的 . 如 人 们 都 知 


结构 ， 许 多 科学 家 认为 蜂窝 的 结构 是 出 于 极 小 原理 : 


m 


道 峰 窜 是 一 种 十 分 类 妙 


这 样 ， 很 
的 六 边 形 


Uses. Jy T DEBER Ode 


尽量 少 ， 就 要 采取 最 经 济 的 设计 方式 . 现在 的 问题 就 是 蜂 帘 是 极 小 曲面 吗 ? 人 们 经 过 多 年 


的 探索 之 后 也 没有 给 出 准确 的 答案 ， 蜂 窝 的 设计 实际 j 
代 的 积累 与 自身 进化 需要 而 建成 的 ， 在 数学 理论 上 很 难 给 日 
造型 ， 从 中 寻求 灵感 来 设计 出 既 美 观 又 实用 的 建筑 . 


人 们 根据 很 多 自然 界 中 的 物理 


要 经 过 时 


上 是 蜜蜂 根据 自己 的 实际 需 


F Plateau 的 贡献 ， 人 


门 习 惯 把 决定 由 空间 Jordan 


| 线 张 成 的 拓扑 类 型 为 


精确 地 揭示 [You,2008] 当然， 


圆 盘 的 极 


小 曲面 问题 称 为 Plateau 问题 [M.Struwe,1988]. 该 问题 早已 成 为 许多 学 者 研究 的 重要 课题 . 


1865 £, H.A.Schwarz 给 出 了 边界 是 边 长 为 1 的 正 立 方 
小 曲面 的 解析 形式 ,这 是 Plateau IET] os — 4 ER 


体 上 的 四 条 边 构 成 扭 四 


边 形 的 极 


解 [Nitsche,1988]. 几乎 同时 ， 另 一 位 德 
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国 数学 家 了 B.Rieman 也 得 到 了 同样 的 结果 [Nitsche,1988]. 1930 年 , J.Douglas fil T. Rado 同 


时 各 自 独 立 的 发 展 了 一 套 解决 Plateau 问题 的 全 新 的 方法 , 他 们 的 方法 很 简便 : Rado 使 用 


共 形 映射 和 逼近 问题 中 的 极限 定理 得 到 Plateau 问题 的 解 [Rado,1930]， 而 Douglas 则 使 用 


变 分 计算 的 直接 方法 [Douglas,1931], 他 也 因此 获得 了 当年 的 第 一 届 菲 尔 兹 奖 .他 们 的 方法 本 
质 上 是 将 面积 极 小 问题 转换 为 Dirichlet 能 量 积 分 问题 , 所 以 Plateau 问题 实质 上 是 决定 给 
定 曲线 的 最 优 参 数 即 等 温 参 数 化 问题 . 参数 曲线 曲面 表示 一 直 是 描述 几何 形状 的 重要 


H, Coons,Bezier 等 对 参数 几何 造型 商定 了 基础 . Bézier 曲线 是 CAGD 中 最 基本 的 造型 工 


具 之 一 .为 了 使 得 造型 能 够 方便 调整 设计 ， 一 般 采 用 通过 调整 控制 顶点 来 调整 造型 的 方法 ， 
这 样 的 调整 比较 简便 . 不 同 的 控制 顶点 就 会 产生 不 同 的 曲面 . 采用 参数 形式 , 一 方面 可 以 突 
出 几何 不 变量 的 含义 ， 另 一 方面 为 过 渡 到 微分 流 形 理论 打 基础 

在 那 之 后 ，Federer 等 人 构造 了 一 条 Jordan 曲线 作为 边界 曲线 ， 在 该 边界 曲线 上 找 不 
到 有 限 亏 格 的 曲面 作为 极 小 曲面 [Fleming,1969,Federer,1969,Federer,1970]. 他 们 的 工作 推 
广 了 Plateau 问题 的 研究 范围 ， 也 说 明了 我 们 不 仅 需要 对 极 小 曲面 的 拓扑 类 型 进行 分 类 ， 


而 且 还 需 借 助 微分 几何 的 概念 以 及 微分 流 形 的 概念 . 更 进一步 , Morrey 后 来 在 Rieman 衬 


间 中 解决 了 更 广义 的 Plateau 问题 [Morrey,1966] 

从 微分 几何 的 观点 我 们 可 以 更 为 精确 地 测量 几何 体 局 部 的 性 质 ， 而 不 是 仅 停留 在 外 部 
形状 的 研究 上 ， 这 使 得 研究 从 现象 到 本 质 产生 飞跃， 使 得 研究 更 具有 科学 性 .为 了 探究 曲面 
的 弯曲 程度 ， 我 们 会 计算 该 曲面 的 绝对 平均 曲率 ， 平 均 绝对 曲率 小 的 曲面 是 较为 平坦 的 ， 
这 里 使 用 曲率 来 刻画 曲面 弯曲 程度 . 平均 曲率 的 概念 有 助 于 进一步 刻画 极 小 曲面 ， 曲 面 在 
某 点 的 平均 曲率 为 0 就 意味 着 曲面 在 该 点 的 两 个 主 法 曲率 大 小 是 相等 的 [Mei, Huang, 2003]. 

极 小 曲面 在 现代 膜 结构 工程 中 有 着 很 重要 的 应 用 . 在 几何 上 ， 根 据 曲 面 在 两 主 曲 率 方 
向 的 曲率 乘积 ,可 将 曲面 分 为 正高 斯 曲面 、 零 高 斯 曲面 和 负 高 斯 曲面 三 大 类 .正高 斯 曲面 的 
两 个 主 曲 率 半 径 均 位 于 曲面 的 同 侧 ， 如 球面 ， 这 类 曲面 也 称 为 同 向 曲面 (synclastic surface). 
零 高 斯 曲面 在 两 个 主 曲 率 方向 中 有 一 个 方向 的 曲率 为 零 , 如 柱 面 . 负 高 斯 曲面 的 两 个 主 曲率 
半径 分 别 位 于 曲面 的 两 侧 ， 如 鞍 面 ， 这 类 曲面 也 称 为 互 反 曲面 (anticlastic surface). 前 两 类 
曲面 较 多 应 用 于 类 似 网 过 结构 那样 的 刚性 结构 中 ， 而 互 反 曲面 则 是 索 网 结构 和 膜 结构 中 常 
用 的 曲面 形式 [Wu,Xu,2004]. 负 高 斯 曲率 曲面 在 建筑 业 、 造 船 业 及 航空 工业 中 都 有 广泛 的 
应 用 .在 建筑 上 ,可 利用 负 高 斯 曲率 曲面 来 设计 马鞍 形 的 房 项 ,如 杭州 市 的 浙江 省 体育 馆 , 其 房 


项 就 是 马鞍 形 曲面 .在 造船 工业 中 ,要 采 
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面 的 研究 是 几何 设计 的 一 个 重 
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以 通过 改 
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La [Jin, Yang,2000]. 


履 盖 较 大 面积 而 不 增加 不 必要 的 重 
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上 一 组 膜 片 来 形成 .设计 几 


变 内 外 支承 的 布置 来 得 到 许多 负 高 
方面 的 
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币 们 往往 根据 极 小 
的 多 方面 


面 的 原型 
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典型 的 负 高 斯 曲面 ， 其 实 可 以 采 


H 
ZNIE 
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计算 机 科学 技术 的 飞速 发 展 ， 使 得 膜 结 


用 


因素 .在 设计 过 程 中 ,支承 构 伯 


E 间 的 跨度 和 承受 的 荷载 相 联系 .这 样 ， 通 过 逐步 设计 最 终 得 到 最 优 的 建筑 


E 来 进行 建筑 设计 ,方案 要 使 
因素 与 其 功能 相 适 应 .最 简单 的 方案 往往 最 有 效 的 .工程 ， 


扬 曙 面 .进行 综合 设计 时 要 同时 考虑 奸 


F 的 布置 和 尺寸 以 及 膜 面 的 形状 应 


其 结构 、 外 观 形态 


党 月 


的 曲面 形式 是 多 种 多 样 的 . 
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用 于 博览 会 、 体 育 场 等 大 型 的 公 


ES 


丹佛 新 


极 小 曲 


构建 筑 更 加 受 
ES E. 如 南宁 澳 海蓝 
际 机 场 候 机 厅 [C.W. Fentress,Jia,1996] 等 . 


面 在 建筑 设计 、 飞 机 轮船 制造 业 、 结 构 化 学 、 材 料 科 学 、 植 物 学 、 动 物 学 等 领 
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中 心 及 美国 


国家 游泳 体育 


域 都 有 着 


E 常 重要 的 应 用 .另外 ， 它 形状 独特 并 


k 有 超凡 的 美感 ， 许 多 艺术 家 收藏 极 小 


作为 模型 


的 作品 ， 


因此 将 极 小 曲面 引入 CAGD 有 着 积极 的 意义 . 


§ 1.2 极 小 曲面 研究 在 CAGD 中 的 发 展 概况 
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边界 的 Bézier 
探索 ， 对 极 小 曲面 的 发 展 做 


AR 


面 中 面积 是 最 小 | 


类 特殊 的 非 线性 PDE 曲面 , 文献 [Man,Wang,2003] 利 月 
# 有 给 定 边 界 的 极 小 曲面 的 B- 样 条 函数 曲 
引入 CAGD 造型 系统 . CAGD 领域 中 的 Plateau-Bézier 问题 可 以 
线 为 边界 的 Bezier 
的 [Jin,1999]. xj 
了 贡献 . 极 小 曲面 是 平均 
圆 偏 微分 方程 的 解 . PDE 曲面 造型 方法 利用 此 类 方程 中 1 


日 非 线性 约束 优 
EMEN, 将 极 小 曲面 正 
述 为 : 给 定 空间 中 的 一 


空间 
面 , 使 得 该 曲面 是 在 所 有 具有 相 
F 这 一 问题 ， 很 多 学 者 进行 了 不 同 角度 
率 为 0 的 曲面 ， 从 本 质 上 是 要 


届 微 分 算 子 实 际 上 是 一 种 


它 部 
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滑 化 过 程 ， 由 这 种 方法 生成 的 曲面 自然 、 光 顺 ， 因 此 被 广泛 应 用 于 船体 、 飞 机 外 形 等 


的 设计 . 


Lluch A.,Monterde,2003] 中 分 别 用 和 矩形 域 上 的 张 量 积 Bézier 


基于 Dirichlet 能 量 函 数 最 小 化 的 方法 ， 在 文献 [Monterde,2004] 及 文献 [Arnal A., 


上 张 量 积 Bézier 曲面 为 调和 曲面 的 充 要 条 件 ， 并且 证 明了 该 张 量 积 
应 的 两 条 边界 曲线 的 控制 顶点 完全 决定 ， 这 给 设计 带 来 了 很 大 的 方便 
利用 满足 等 温 参 数 条 件 的 调和 参数 曲面 是 极 小 曲面 这 一 性 质 ， 很 多 学 者 对 给 定 边界 的 


[Monterde,2003] 4] [Monterde,2004] F 25 JE FA i A 


面 和 三 角 域 上 的 B-B 曲 
以 表示 极 小 曲面 .但 是 这 里 的 方法 都 很 复杂 , 在 实际 应 用 中 极 不 方便 .实际 上 , 如 果 一 张 参数 
面 满足 等 温 参数 条 件 ， 那 么 它 为 极 小 曲面 的 充 要 条 件 是 它 为 调和 
I 曲面 来 近似 表示 极 小 曲面 ， 给 出 了 和 


面 来 近 


pun 


极 小 曲面 的 造型 设计 做 了 研究 . Farin 和 Hansford [Fairn,Hansford,1999] 认 为 直接 用 平均 


为 0 来 刻画 Plateau-Bézier 问题 ， 在 造型 设计 时 存在 很 多 不 便 . 他 们 ] 
提出 了 离散 拉 普 拉 斯 算 子 生成 网 格 方法 ， 这 种 方法 计算 量 少 ， 生 成 的 曲面 具有 保 形 性 
量 泛 函 最 小 原则 . 虽然 他 们 的 这 种 方 济 


直接 针对 控 人 


面 .文献 
E 形 域 
面 的 控制 网 格 由 其 相 


曲率 


草 顶 点 本 身 


及 能 


显得 过 于 粗糙 ， 但 是 他 们 也 开创 了 研究 


Plateau-Bézier 问题 的 新 方法 . 在 那 之 后 ，Monterde [Monterde,2003] 深 入 研究 Dirichlet 函数 
的 性 质 ， 研 究 了 特殊 的 一 类 Plateau-Bézier 问题 : 其 给 定 的 边界 曲线 是 Bézier 曲线 .他 通过 
用 Dirichlet 泛 函 来 奉 代 面 积 泛 函 来 通 近 Plateau-Bézier 问题 的 解 .接着 Monterde 和 Ugail 


[Monterde,Ugail2004] 给 出 了 基于 椭圆 形 侦 微 分 算 子 的 给 定 边界 来 构造 曲面 的 新 方法 , 特别 


和 情形 的 一 般 化 . 


究 了 用 四 阶 PDE 的 方法 来 产生 给 定 控 


地 ， 致 力 于 调和 及 双 调 和 Bézier 曲面 的 设计 ， 从 中 得 出 了 重要 结论 : 任何 双 调 和 的 B 
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由 边界 的 Bézier 曲面 .这 种 方法 是 对 前 面 调和 及 双 调 


ézier 


与 此 同时 ， 满 家 巨 和 汪 国 昭 研究 了 B 样 条 极 小 曲面 造型 [Man,Wang,2003], 他 们 利用 等 
温 参数 极 小 曲面 的 调和 性 质 , 证 明了 三 次 等 温 参 数 多 项 式 极 小 曲面 的 唯一 性 .另外 ， 徐 岗 和 


汪 国 昭 利用 方向 导数 研究 了 三 角 域 上 的 调和 B-B 曲面 [Xu,Wang,2006]， 给 出 了 三 角 域 


上 的 


B-B 曲面 为 调和 曲面 的 充 要 条 件 ， 并 且 证 明了 任何 一 个 三 角 域 上 的 调和 B-B 曲面 的 控制 网 


格 均 由 它 的 第 一 层 和 第 二 层 控制 顶点 决定 ， 最 后 还 对 极 小 曲面 在 建筑 设计 中 的 应 用 进行 


初步 探讨 . 这 些 论 文 的 研究 探索 对 进 


步 研 究 极 小 曲面 的 造型 具有 习 


尽管 极 小 曲面 造型 设计 在 计算 机 辅 


助 几 何 设 计 领 域 取得 了 很 多 重要 成 果 ， 但 是 其 


EE 要 作用 . 


行 了 


研究 


的 范围 还 是 比较 狭窄 ， 并 不 能 完全 满足 实际 的 需要 . 我 们 知道 Bézier 曲线 不 能 精确 表 
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多面 外 的 二 次 曲面 ， 而 


里 调和 Bézier 


浙江 大 学 硕士 学 位 论文 
抛物 线 以 外 的 圆锥 曲线 ， 而 有 理 Bézier 曲线 既 能 表示 多 项 式 曲 线 (如 飞机 机 身上 的 纵向 曲 
线 )， 又 能 表示 圆锥 曲线 (如 飞机 机 身上 的 横 截 面 线 )， 可 以 把 两 者 统一 处 理 ，Bkzier 
有 理 Bézier 曲线 的 特殊 形式 . 类 似 地 , Bézier 曲面 不 能 精确 表示 除 抛 4 
AH Bézier 曲面 能 克服 这 一 缺陷 . 因此 把 以 往 构造 极 小 曲面 的 方法 推广 到 有 
面 的 设计 有 着 很 广泛 的 应 用 价值 . 
$13 本 文 的 主要 工作 


基于 以 上 原 


有 限 差 分 法 来 构造 出 


骨 和 曲面 . 基于 有 


如 果 直接 利 


Bézier 曲面 ， 处 理 过 程 也 是 很 复杂 的 ， 
程 的 形式 , 本 文采 用 了 有 


限 大 时 ， 


j 调 和 算 子 进行 求解 ， 则 会 产生 


再 将 求解 出 的 点 连接 起 来 便 可 产生 所 要 求 的 


| 这 些 点 所 生成 的 网 格 面 就 无 限 接近 于 所 求 的 有 到 


系列 的 高 度 非 线 怕 


因 ， 本 文 主要 探索 了 有 理 调和 Bézier 曲面 的 构造 方法 ， 其 主 


要 


4p 


思 


H 


想 是 


ayy) 


采用 


里 Bézier 曲面 的 复杂 形式 , 求 导 之 后 形式 更 为 复杂 ， 
的 方程 ， 即 使 是 二 阶 有 理 


这 点 可 以 详 阅 后 记 . 为 了 将 解决 方案 转化 为 线性 方 


用 于 任意 nn 阶 有 理 Bézier 情形 ， 并 且 求 解 过 程 都 是 解 线性 方程 组 


A 


本 文 是 这 样 来 组 织 的 : 第 


际 工程 


te E a a 


Ha on 


是 对 本 文 所 需要 的 微分 几何 的 知 


HET 
中 的 应 用 , 其 次 是 极 小 曲面 在 CAGD 领域 中 的 发 展 概况 , 最 后 综述 本 文 的 主要 工作 


用 的 造型 理论 加 以 简单 的 回顾 ， 


为 之 后 的 探索 做 知识 铺垫 .第 三 章 是 有 理 调和 Bézier 曲面 造型 


了 国 


限 差分 的 方法 , 将 要 构造 的 曲面 加 以 离散 从 而 得 出 曲面 上 的 点 值 ， 


面 . 另外 ， 这 是 一 种 逼近 ， 当 划分 的 份 数 无 


调和 Bézier 曲面 . 这 种 方法 适 


际 上 对 于 极 小 曲面 的 研究 及 极 小 曲面 在 实 


识 及 调和 Bézier 


i 


行 的 有 限 差 分 法 .多 


长 


n" 


四 章 主 要 是 应 用 大 量 的 实例 来 验 订 


的 总 结 和 对 将 来 工作 的 展望 . 


FE 这 种 造型 方 没 


的 探索 ,而 


.第 五 章 是 对 


究 了 针对 曲面 来 进 


文 内 容 
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$ 2.1 极 小 曲面 与 调和 曲面 的 关系 


民 多 学 者 致力 于 研究 Plateau-Bézier 问题 , 其 基本 的 解决 方法 都 是 联系 调和 曲面 来 进行 
究 .那么 极 小 曲面 与 调和 曲面 之 间 到 底 存 在 什么 关系 呢 ? 这 一 节 将 对 此 加 以 阐述 . 


eS 


zn 


假定 空间 尺 中 的 参数 曲面 表示 如 下 : 


r(u,v) =(x(u,v), y(u, v), Z(U,V)),U € (—%, +00), v € (—00, +00), 
它 的 第 一 基本 形式 为 


E =<r r, >, F =<r,,r, >,G=<r,,r,>. 


Eer lur, Aru, v) 的 一 阶 偏 导 数 ，<,> 为 向 量 的 点 积 运算 算 子 . 


它 的 第 二 基本 形式 为 


了 =( 人 7 -(rrr»N-(r,r,). 


Hor, sr, Mr, Arlu, v) 的 二 阶 偏 导数 ，(,,) 为 向 量 的 复合 运算 算 子 .曲面 (wy) 的 平均 


曲率 H 和 高 斯 曲率 天 分 别 为 


_EN-2FM+LG , LN-M* 
2EG-F*) ' EG-F?- 


定义 2.1.1 d ilf (u,v) iie E = G, F — 0, WEK HH r (u, v) 为 等 温 参 数 曲 面 . 


定义 2.1.2 Er HH IRE r (us v) Er, +r, = 0, 则 称 曲 面 r(u,v) 为 调和 曲面 . 


定义 2.1.3 FH ru, v) IAE XS iss H HOUSE, WERE IR rn (us v) 为 极 小 曲面 . 


引 理 2.1.1. 等 温 参 数 曲 面 为 极 小 曲面 的 充 要 条 件 是 它 为 调和 曲面 . 


证 明 : 假 设 曲面 r(u,v) ASS we hh, BE =G,F =0, 其 中 


E =<r (u,v),r (u,v) >, F =<r,(u,v),r,(u,v) >,G =< r (u,v),r (u,v) >. 


对 上 述 两 式 两 边 同 时 求 导 可 得 
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<r (u,v), r (U,V) »-«r, (u,v).r (U,V) >=— <r, (u,v),r, (u,v) >. 
即 有 


<r (UV),T,, (U,V) tr, (u,v) >=0 


<r (u,v),r (UV) tr, (u,v) >= 0 


成 立 . 于 是 推 知 向 量 r,, Qu v) +r, (u,v) 平 行 于 曲面 的 单位 法 向 量 N ; 


同时 ， 以 e, f, g 表示 正则 曲面 的 第 二 基本 形式 的 系数 ， 则 有 


Hil-2/C+tsE_letsg 
2 EG-F? 2 E` 


因此 有 


H =0&e+g =0&<N,r (u, v)> +< N,r (u, v)>=0 < N,r tr,»-0 


> 


= = ur EG 
rtr, | =O0or, +r, =0. 证 平 


由 于 直接 求解 给 定 空间 闭 曲 线 上 的 极 小 曲面 是 相当 困难 的 ， 直 到 如 今 ， 也 没有 很 好 的 
形式 解 . 这 时 ,我 们 就 要 转变 思路 寻找 其 他 的 方法 来 求 取 极 小 曲面 . 基于 等 温 参 数 曲面 是 极 


小 曲面 的 充 要 条 件 是 其 是 调和 的 这 一 结论 ， 对 于 Plateau-Bézier 问题 ， 利 用 Bézier 曲面 的 
特殊 性 质 ， 可 以 得 到 很 好 的 极 小 曲面 造型 算法 ， 这 点 在 以 后 的 章节 中 会 加 以 介绍 . 这 种 造 


型 方法 在 理论 上 易 懂 ， 实 现 起 来 也 很 简单 ， 对 极 小 曲面 的 研究 有 重要 影响 . 


IN 


$ 22 Bézier 调和 曲面 设计 


基于 上 节理 论 知 识 ， 这 时 Bézier-Plateau 问题 可 以 这 样 描述 : Zaz BIS Fi lees 


c— 


边界 为 Bezier 曲线 ), 目标 是 要 寻找 基于 给 定 边界 的 面积 最 小 的 曲面 . 求 极 小 曲面 的 问题 就 
可 以 转化 为 求解 能 量 泛 函 的 极 值 . 


Bog C: 空间 中 一 条 闭 的 可 求 长 的 简单 Jordan 曲线 ，Q :曲线 CC 在 平面 上 的 投影 所 围 
成 的 区 域 . 定义 函数 空间 : 


X: {z=2(x,y):Q>R| z(AQ)=C}. 


Plateau 问题 就 是 在 在 中 寻找 面积 最 小 的 曲面 ， 而 曲面 的 面积 可 以 表示 为 
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A(z) = ff iz +z; dxdy. (2.21) 
Q 


XX ^ f EZ KAE A) TE, EAL Euler — Lagrange 方程 : 


o — & E o ay =0 (2.2.2) 
Ox Jitz +z oy Jitz +z, l 
上 述 方程 (2.2.2) 被 称 为 极 小 曲面 方程 . 该 方程 是 非 线 性 的 偏 微 分 方程 , 这 一 问题 的 求解 是 很 
困难 的 .在 探究 Plateau — Bézier 问题 时 ， Monterde [ Monterde ,2003; Monterde ,2004] 


H Dirichlet 的 方法 进行 逼近 来 构造 出 极 小 曲面 . 


% 


II 


[E E E USAT r(u v) :U — R^, IHU =[0,1]x[0,1], 则 其 面积 可 以 表示 为 


r, xr, 


1 
A(r) = | Hudv = | (EG - F^) dudv. (2.2.3) 
U U 


另 一 方面 Dirichlet 泛 函 可 以 表示 为 : 


2 


T+, 


pr)=|( 


^Mudv = | (E(u,v) + G(u,v))dudv. (2.2.4) 
U 


于 已 =<7 (u,v),r, (u,v) >, F =<r,(u,v),r,(u,v) >,G =<r (u,v),r, (u,v) >, 


E20,F 20,G20. 
成 立 . 则 


1 1 
(EG - F?y <(EG)? < EE 


所 以 对 于 任意 的 控制 网 格 及 其 曲面 r(u,v) EA Ar) € Dr) RAHE =G, F =0 


时 等 号 成 立 .而 由 前 面 知 识 可 知 ， 满 足 E=G,F =0 的 曲面 是 等 温 参 数 曲 面 .对 于 等 温 参 数 


曲面 ， 使 得 所 构造 的 曲面 面积 泛 函 A(z) = (f V1+z? +z?dxdy 最 小 等 价 于 使 得 Dirichlet 
Q 
i 


2 
十 


r, 


r, 


D(r)= K ^Mudv = TE t G(u,v))dudv 
U U 
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最 小 .由 微分 方程 的 知识 可 知 ， 要 使 Dirichlet 泛 函 最 小 需 满足 Euler — Lagrange 方程 : 


O’r(u,v) 7 O’r(u,v) 


D AP Ar =0. (2.2.5) 


我 们 将 Dirichlet 泛 函 取 得 极 值 所 满足 的 Euler — Lagrange 方程 定义 为 拉 普 拉 斯 算 


d. 类 似 于 调和 曲面 的 定义 ， 也 可 以 利用 Euler — Lagrange 函数 定义 标准 的 双 调 和 算 子 . 


那么 如 何 构造 给 定 的 控制 边界 是 Bézier 曲线 的 极 小 Bézier 曲面 呢 ? 


首先 ， 假 设 以 空间 中 的 点 瑟 为 控制 顶点 的 Bézier 曲面 表示 为 : 


R(u,v) = Syg (u) Bi" ()P;, HF R :[0,1]x[0,1] > R°. 


i=0 j=0 


定理 2.2.1 [Monterde, Ugail, 2004] 4s R? HRI, {P} 。， 它 所 对 应 的 Bezier 


Wi R:[0,1]x[0,1] > R? Æ M hm, 即 ARU, v) =0 的 充 要 条 件 是 对 于 任意 的 


iefl, nn}, je {l,m} 有 以 下 方程 组 成 立 . 


Pa nio tP uua 117 24,49) 

+P iy (ap i-11 72an i-2,2) + P,, ja ni-22 
+P i, +24, j,0 +P ija am j- = 2a, ;0) 
+P, ja (a 
+P (a, io = 24, iy * 0,153 
2a 


eg AA 24, ;22) +P i, j-24 m, j-2,2 


TO. 5,0 一 Q, jaa + Gn, j- 22) = 0. 


其 中 ie {0,…,n 一 2}， 
-(n-i)(n-i-1),a,, -2(--1D)(n-i-1),a,, =(i+D(i+2). 
"igll-snnjell-um] It, P.-0; 
"Aig[0,-.n-2]H., duo = Gy = au = ORE. 


证 明 : 


将 拉 普 拉 斯 算 子 作用 于 Bézier tiii R(u,v) 可 得 : 
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o 0 
AR(u,v) = E 十 Z jrun 


n-2 m n- m 
-n(n-15 Br WBF WAP, 


n m-2 n vies 
mam -1)» , 25 B; (u) B? ^ (v)A"? P, 


ij? 
APP -pP .—2P +P 
GÍ = i,j+l DI^ 


i,j42 


i-2,j 


2,0 x 
A P, = 五 —2P.., tP, 


N 
4| 


为 了 简化 上 述 多 项 式 ， 我 们 需要 引入 以 下 的 变换 : 
(n-iX(n-i-DB; (1) 


(t) - 6 - DG--2)B; .(t)). 


n-2 A 1 
RESINE n(n—1) 


32i D) (ni -1) B2, 


b, Cp 分 别 定义 为 ; 


对 于 ie {0,…,n 一 2} , 我 们 把 a;， 
a, = (n-i(n—i-1),b,, = 2 +1)(n-i-1),c¢,, = G+ DG+2), 


x Tie(o-sn-2jM, a, —b, = cn =0. 


则 有 


AR (u,v) = > 2 ,Br (Bi) 
(n(n m D(a,A™P, 2 Dey Ps m Ga ANP 


+m(m—INa mA” Pj +b, , AV P ire LAT P a). 


面 ， 不 妨 记 其 控制 顶点 为 


面 ， 仍 然 是 一 张 Bézier 


经 拉 普 拉 斯 算 子 作用 过 的 Bézier 
ETTAN: AR(u,v) 20e {2 4 =0. 


lO ， 则 利用 Bézier 基 的 线性 无 关 性 


将 上 述 得 出 的 等 式 加 以 整理 化 简便 可 以 证 明 此 结论 . 


当 m =n 时 ， 定 理 2.2.1 [Monterde, Ugail, 2004] 的 结论 简化 为 以 下 结果 
(b; ,, E 2a,,) is P ba, Nd 2€; 54) 
i424 jm +P, a (bias X. 2a im) (2.2.7) 


+P ja (Diin g 2€, 1m) T P, -2C s 


+P, (a, — 2b in 6; 4, + jm — 2b, 


0-P., a, +P 


i+2,j “in 


T ce RP 


i-2,j “i-2,n 


m F C) 


由 顶点 之 间 有 着 好 的 性 质 ， 即 其 


道 ， 对 于 Bézier 调和 曲面 ， 控 人 
Hl) JU p 


由 上 述 定理 我 们 可 以 知道 
中 的 一 些 控制 顶点 是 其 他 一 些 控制 顶点 的 线性 组 合 .精确 的 说 ,第 一 排 和 最 后 一 排 控 人 
可 以 完全 确定 整个 Bezier 调和 时 面 .这 一 点 是 如 何 证 明 的 呢 ? 首先 看 一 下 下 面 的 引 理 
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引 理 2.2.1 [Monterde,2003] 3 f(u.v)= $^. , auu v 是 n 次 调和 多 项 式 ， 其 中 


n>2， 则 有 以 下 结论 成 立 : 


(D 如 果 n 是 奇数 ， 则 
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所 有 的 系数 人 ay} _， TAH faya} 完全 确定 ， 即 所 得 


k=2,1=0 


系数 都 可 以 表示 为 第 


(2) 如 果 nn 是 偶数 ， 则 所 有 的 系数 {ay} 


定 ， 即 lau], 


k=2,1=0 


行 及 第 二 行 系数 的 线性 组 合 . 


和 au 可 以 | 


k=2,1=0 


Ma, PUVA {ay}, (ag), 线性 表示 


1=0 


我 们 知道 Bézier HES HSE [8] Ze wp ULRICH, 这 两 组 基 


以 得 到 以 下 的 推论 . 


推论 2.2.1 [Monterde, 


一 一 


o, bo 及 {a ju 完全 确 


a 


Ugail, 2004] 18% X (u, v) = Dae 


和 的 Bézier 曲面 ， 其 控制 项 点 {P AUC FA 


HERR, BNP, fe 


k,l=0 


是 等 价 的 . 由 引 理 2.2.1 可 


Bi (u) Bj (v) B, Æ n Ér il 


-LI-0 可 以 | Usa EU 完全 


f 


O 如 果 n 是 奇数 , 中 间 的 n 一 1 排 控制 顶点 可 以 由 第 一 行 和 最 后 一 行 的 控制 顶点 线 


zi 


dE. 


(2) 如 果 nn 是 偶数 ， 中 间 的 控制 顶点 {PB} 和 角 点 ,可 以 由 第 一 行 和 最 后 一 


FTR {Pu} UP, 


有 了 以 上 的 理论 知识 做 准 


小 一 完全 确定 . 


备 ， 我 们 现在 要 解决 的 问题 是 : 


的 边界 ， 已 知 两 条 相对 边 
面 呢 ? 


给 定 n Bézier 曲线 组 成 


N 


i=0 j=0 


即 二 阶 情形 ， 这 时 已 知 空间 中 的 控制 顶点 矩阵 { 己 Si pa 


民 据 定理 22.1 的 特殊 情形 (m 


P P P» 


00 01 


* 


* * 
Py P, P, 


界 的 控制 顶点 ， 那 么 如 何 求 得 在 该 给 定 边 


首先 , 看 一 下 简单 情形 , 对 于 曲面 R(u,v) = Y Y B"0DB"O)P 


界 上 的 调和 Bézier HH 


当 n = 2,m = 2 时 ， 


ij? 


=n) 所 得 到 的 等 式 可 得 到 以 下 等 式 : 
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1 
Py, xp 2 Q. Py 十 Py 一 2P, 十 P). 


1 
ju 00 + Py + Pao P5), 
i (2.2.8) 
Py, = a 2 (Po 未 2P, ES 2P, + P), 
1 
P, ic xi 00 + Py +2Po -Pao P 2). 
Be) Hy WAAAY A (2.2.8) P BY DG AZ PEZ; Be OR th PG ARAB Po Pi P.. t 


而 可 以 构造 出 所 要 求解 的 调和 曲面 . 


我 们 再 来 看 下 简单 的 三 次 情形 : n = 3,m = 3 时 . 


这 时 已 知 室 间 中 的 控制 顶点 矩阵 {P, bP OM 
Po Py, Py Py 


* * * * 
* * * * 


pop P. 


这 时 候 上 只 知道 第 一 排 和 最 后 一 排 控制 项 点， 这 就 需要 一 个 含有 8 个 方程 的 方程 组 ， 从 中 反 


求 出 未 知 的 顶点 . 根据 定理 2.2.1 的 特殊 情形 (m = n) 所 得 到 的 等 式 组， 我 们 可 以 得 到 以 下 


等 式 组 : 


1 
P, = a 十 2803 +2P, + P.) 


1l 
ig I i 
1 

=5 Pn + 4 + Po +2P,), 


P. E +2P +2P,, +4P,,), 


30 


1 (2.2.9) 
Po = 3 z CABG - AP, +2Py +2P — 2P,, + P 3), 


1 
P = 3 270 十 Py, + AP. — AP, +2P 32) 


1 
P, = Bo 4P, T AP, ns P. -2P 3+ 2P,;), 


1 
Py, = 3 (n —2P oy + 2Po; +2P, - AP 3+ 4P,3). 


从 (2.2.9) 中 我 们 可 以 反 求 出 未 知 控制 项 点 ， 且 其 解 是 唯一 的 . DKA Ae nx BERBA 
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呢 ? 答案 是 肯定 的 , 我 们 可 以 从 定理 2.2.1 中 得 到 求解 的 算法 . 即 给 定 控制 边界 的 第 一 排 和 


最 后 一 排 控制 项 点 ,我 们 便 可 以 通过 曲面 调和 这 一 条 件 来 确定 出 调和 
那么 它 便 是 极 小 曲面 . 这 一 算法 对 于 


如 果 该 曲 


造型 设计 


^m 


FY AS AMT BA 


方程 组 ， 将 得 到 在 每 个 结 点 处 近似 于 真 值 的 函数 值 ， 然 后 再 利用 所 


得 到 定 解 


区 域 作 网 
为 了 保证 


面 是 等 温 参数 曲面 ， 
有 重要 的 意义 


差分 方法 是 用 于 微分 方程 定 解 问题 求解 的 最 
差分 法 的 基本 概念 是 用 差 商 代替 微 商 ， 从 而 把 1 


§2.3 有 限 差 分 法 介绍 


面 所 有 的 控 种 


AES 


J Bézier 调和 


局 微分 方程 转换 为 差分 方程 . 将 所 研究 


常用 的 数值 方法 [Yu, Yang,2003]. 


IN 


数 的 域 细 分 为 许多 网 格 ， 然 后 列 出 对 应 每 个 交叉 点 ( 结 点 ) 的 差分 方程 ， 解 该 差分 


问题 在 整个 区 域 上 的 近似 解 [An,Liu,Guo,1996]. 


有 限 差分 法 的 主要 内 容 包 括 [http://baike.baidu.com/view]: 如 何 


§ 值 方法 便 可 以 从 离散 解 


民 据 问题 的 特点 将 定 解 


格 剖 分 ;如何 把 原 微分 方程 离散 化 为 差分 方程 组 以 及 如 何 解 此 代数 方程 组 .此 外 
计算 过 程 的 可 行 和 计算 结果 的 正确 ， 还 需 从 理论 上 分 析 差 分 方程 组 的 性 态 ， 包 括 


解 的 唯一 性 、 存 在 性 和 差分 格式 的 相 容 性 、 收 敛 性 和 稳定 性 .对 于 一 个 微分 方程 建立 的 各 种 


差分 格式 


要 求 . 另外 ， 
的 解 ,这 就 是 收敛 性 


， 为 了 有 实用 意义 ， 


值 ， 直 到 与 初始 值 有 


HEE n 层 的 近似 
门 到 后 面 各 层 的 值 ， 如 果 误 将 的 


个 基本 要 求 是 它们 能 够 任意 通 近 微分 方程 ， 这 就 是 相 容 性 
个 差分 格式 是 否 有 用 ， 最 终 要 看 差分 方程 的 精确 解 能 否 任意 逼近 微分 方程 
的 概念 . 此 外 , 还 有 一 个 重要 的 概念 必须 考虑 , 即 差分 格式 的 稳定 性 因 
为 差分 格式 的 计算 过 程 是 逐 层 推 进 的 ， 在 计算 第 n +1 层 的 近似 值 时 要 月 
关 . 前 面 各 层 若 有 售 入 误差 ， 必 然 影 ! 


影响 越 来 越 大 ， 以 致 差分 格式 的 精确 解 的 面貌 完全 被 拖 盖 ， 这 种 格式 是 不 稳定 的 ， 相 反 如 


果 误 差 的 


传播 是 可 以 控制 的 ， 就 认为 格式 是 稳定 的 .只 有 在 这 种 情形 ， 差 分 格式 在 实际 计算 


中 的 近似 解 才 可 能 任意 通 近 差分 方程 的 精确 解 . 


关于 差分 格式 的 构造 一 般 有 以 下 3 种 方法 [http://baike.baidu.com/view]. 最 常用 的 方法 
分 法 ， 比 如 用 差 商 代替 微 商 等 . 另 一 方法 叫 积分 插值 法 ， 因 
疗程 常常 反映 物理 上 的 某 种 守恒 原理 ， 一 般 可 以 通过 积分 形式 来 表示 . 此 外 还 可 以 用 


是 数值 微 


为 在 实际 问题 中 得 出 的 
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这 种 方法 在 解决 plateau — Bézier 问题 中 有 着 广泛 的 应 用 , 基于 Bézier 调和 曲面 的 特 


性 来 求解 出 相应 的 调和 曲面 ，Farin 1 Hansford 4 (Farin, Hansford 1999] 中 应 用 这 一 


方法 ， 有 具体 描述 如 下 : 


Ben vk Bézier SHH THR RW R(u,v) = 》》,B”(u)B”(v)P, 要 构造 出 极 小 曲 


i=0 j=0 


面 ， 需 要 使 得 能 量 泛 函 最 小 .特别 地 对 于 Coons 曲面 片 ， 我 们 知道 其 能 量 泛 函 为 
| R? dudv » JUPU =[0,1]x[0,1] 
U 


要 使 其 能 量 泛 函 最 小 ， 需 使 其 满足 Euler- Lagrange PDE: 


=0. 


uuvv 


如 果 我 们 将 此 极 小 原理 应 用 到 离散 的 Coons 曲面 片 ， 我 们 可 以 得 到 以 下 结论 : 


m-l n-l 


i=0 j=0 


Jeh ANB, =b b, —b, +b.. 


i+l j+ “iH, j vij ij 


下 面 我 们 来 看 一 下 有 限 差分 方法 的 简单 应 用 . 


我 们 将 拉 普 拉 斯 算 子 应 用 到 给 定 的 泛 函 六 [u,v] > R ,在 有 限 差分 网 格 上 用 i,j 来 定 


义 网 格 节 点 ， 该 网 格 步 长 为 hh ， 则 节点 的 坐标 可 以 表示 为 : u=ih,v= 六， 可 以 得 到 


> > > > — 一 


eh = ie 十 at = EH S 0 (2.3.1) 
从 (2.3.1) 等 式 中 我 们 容易 知道 
一 1 > > > > 
Xi,j peu Xi, jt Xi, jud "a (2.3.2) 
从 (2.3.2) 式 中 可 以 看 出 中 间 的 未 知 控制 顶点 可 用 其 周围 的 控制 顶点 来 表示 ， 那 么 是 否 存 在 
特殊 的 模版 来 解决 这 一 问题 ? 下 面 我 们 将 引入 mask 定义 . 
如 果 以 调和 mask 的 形式 来 表示 (2.3.2) 为 
0 1 0 
—] e i] (2.3.3) 
0 1 0 
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这 里 mask 被 定义 为 是 用 来 计算 中 间 未 知 点 的 模板 . 从 几何 上 可 以 这 样 说 ， 未 知 控制 顶点 
可 以 由 其 上 下 左右 四 个 控制 顶点 来 确定 ， 这 样 就 大 大 减少 了 计算 量 . 


三 章 有 理 Bézier 极 小 曲面 设计 


$3.1 AH Bézier 极 小 曲面 造型 问题 的 提出 
有 理 Bézier 曲线 在 几何 造型 中 有 其 独特 的 优势 ， 一 方面 ， 可 以 精确 表示 圆锥 曲线 ， 另 
一 方面 ， 引 入 了 权 因 子 ， 使 得 几何 造型 的 调整 更 加 灵活 . 类 似 地 , Bézier 曲面 不 能 精确 表示 
除 抛物 面 外 的 二 次 曲面 , 而 有 理 Bézier 曲面 能 克服 这 一 缺陷 所 以 ， 将 Plateau-Bezier 问题 
推广 到 有 理 形式 是 必要 的 . 


ig (P,)77 为 加 空间 中 (n+Dx(m+DD 个 点 向 量 , 4 wa 27 为 权 因 子 集 ， 则 由 此 生 


k,l=0 k,l=0 


RET] n x m 次 有 理 Bézier 曲 面 可 表 为 : 


Bi (u) Br (v) wy Pu 
R = R(u,v) - E99 — — — — —,  uveU =[0,1] [0,1]. (3.1.1) 
B, (u)Bj (v) wy 


k=0 1=0 


ii R(u,v) 的 第 一 基本 形式 之 系数 为 E,F,G 则 与 此 曲面 相对 应 的 面积 函数 [do 


Carmo,1976] 定 义 为 : 


A(R) = | VEG — F’ dudv. (3.1.2) 
U 


通过 (3.1.1) 可 知 有 理 Bézier 曲 面 是 由 其 控制 顶点 及 其 权 因子 唯一 确定 的 . 因此 可 以 这 样 来 描 
述 Plateau-Rational-Bézier 问 题 :给 定 空间 中 一 条 闭 的 边界 ， 在 该 边界 上 寻找 面积 最 小 的 有 理 


Bekzier 极 小 曲面 .这 一 问题 可 以 等 价 的 描述 为 :给 定 边 界 的 控制 顶点 和 权 因 子 ， 在 所 有 具有 上 


述 相 同 边界 控制 顶点 和 权 因 子 的 有 理 Bkzier 曲 面 中 ,寻找 面积 泛 函 ACR) 达到 最 小 的 曲面 . 


A(R) 中 涉及 到 参数 曲面 的 第 一 基本 形式 ,了 『,G ， 很 显然 需要 对 (3.1.1) 进 行 求 导 . 由 于 
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有 理 多 项 式 求 导 是 十 分 复杂 的 ,求解 过 程 涉及 到 高 度 非 线性 方程 求解 . 这 时 应 该 转换 思路 从 


其 它 角度 来 解决 这 一 问题 . 能 和 否 通 过 求解 线性 方程 组 来 解决 呢 ? 


从 第 二 章 中 可 知 ， 当 EE= G,F =0， 考 虑 曲面 的 Dirichlet 函数 


(E+G) 
U 


及 如 何 使 得 (3.1.3) 取 极 值 ， 当 A(R) = DCR) 成 立时 所 求 得 的 曲面 就 是 极 小 曲面 . 


值得 注意 的 是 : 曲面 的 Dirichlet 函数 依赖 于 曲面 的 参数 化 形式 , 它 的 值 会 随 着 不 同 的 
参数 化 形式 而 改变 . 但 是 ， 曲 面 的 面积 函数 不 依赖 于 曲面 的 参数 化 形式 , 它 是 曲面 的 内 部 
几何 量 . 因此 ,前 者 只 能 作为 后 者 的 一 个 近似 度量 . 换 名 话说 ,在 曲面 的 等 温 参 数 化 形式 下 ， 


Dirichlet 函数 的 最 小 值 等 于 曲面 面积 函数 的 最 小 值 . 


由 偏 微分 的 知识 可 知 : 要 使 (3.1.3) 式 取得 极 值 ， 应 使 其 满足 Euler — Lagrange 方程 : 


R,,,+8,, =9. 


这 就 是 我 们 经 常见 到 的 调和 方程 ,其 求解 模式 比较 固定 , 这 一 理论 知识 在 极 小 曲面 的 造型 过 


程 中 有 着 很 广泛 的 应 用 . 


§ 3.2 有 理 Bézier 极 小 曲面 造型 理论 知识 


对 于 一 张 n xn 有理 Bézier 曲面 : 


22, Br 0B’)w,P, 
R=" (3.2.1) 
> >》B (WB’ (yw, 
i-0 j=0 
为 了 表示 方便 ， 记 其 分 子 ， 分 母 分 别 为 
X =9_ BWB’ vw,P, 
5 (3.22) 
7 = > >》B (WB’ (vw, 
i=0 j=0 


要 构造 出 有 理 Bézier 调和 曲面 ， 应 先 对 有 理 Bézier 曲面 进行 求 导 : 
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X Y-Y X 
R, = u y? u 
R= DEDE 


(X 了 -了 X)r -2Y Y Y +2XYY? 


R, = Y^ 
pa 4 -Y X)Y? -2Y Y, X, - 2XYY; 
vw Y? : 


从 上 节 内 容 可 知 , 要 使 得 其 面积 泛 函 ACR) 最 小 , 则 曲面 应 满足 Euler — Lagrange 77 


程 uu +R, = 0, 即 


(Xa + X, )Y? - XY(Y,, «Y, ) -2Y (X,Y, + X,Y,)+2X(¥7 +Y?) 
Y? 


= 0. 


(X, + XY -XYY „+Y „)-2Y(X, Y, +X Y)422X (Y) &Y2)) 20. (323) 


对 于 任意 的 调和 的 有 理 Bézier 曲面 都 必须 满足 (3.2.3) 式 ， 现 在 的 问题 就 是 如 何 根据 
(3.2.3) 来 求解 呢 ? 也 就 是 如 何 用 已 知 的 边界 控制 顶点 和 权 因 子 来 表示 出 整个 有 理 Bézier 
面 呢 ? 

当然 ， 最 直接 的 方法 是 根据 (3.2.3) 对 (3.2.1) 直 接 进 行 求 导 ， 利 用 Bézier 基 的 线性 无 关 
性 来 求 出 控制 顶点 之 间 的 关系 . 通过 探索 ， 发 现 这 一 过 程 是 很 困难 的 ， 问 题 在 于 有 理 多 项 
式 求 导 后 次 数 会 增加 ， 通 过 升 阶 利用 基 的 线性 无 关 性 来 求解 就 很 困难 . 


观察 (3.2.3) 发 现 , 其 中 有 一 项 是 关于 分 母 的 , 如 果 分 母 是 调和 的 , 那么 也, + 了 ,,)=0， 


(3.2.3) 就 简化 为 


(Xa +X „Y? -27(X,Y, (XY FONG ey) =: (3.2.4) 


另 一 方面 ， 欲 求 曲面 的 分 母 相当 于 多 项 式 Bézier 曲面 的 每 一 个 分 量 ,其 调和 模式 可 以 轻松 
地 求 出 . 这 是 因为 对 于 调和 Bézier 曲面 ， 已 经 有 很 好 的 性 质 和 结果 ， 如 给 定 控制 顶点 的 第 
一 排 和 最 后 一 排 ， 其 他 控制 顶点 可 以 由 第 一 排 和 最 后 一 排 的 控制 项 点 来 线性 表示 ， 从 而 可 
以 构造 出 整个 调和 的 Bézier 曲面 . 

那么 可 以 采用 分 而 治之 的 方法 吗 ? 先 使 分 母 调和 ， 然 后 再 达到 整个 有 理 Bézier 曲面 调 


和 . 这 是 全 部 有 理 Bézier 调和 曲面 中 的 重要 一 类 . 现在 的 问题 是 求 出 的 曲面 是 否 比 一 般 的 
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曲面 性 质 更 好 呢 ? 如 绝对 平均 曲率 较 小 ， 更 加 逼近 于 有 下 
因子 的 选取 . 对 于 有 理 Bézier 曲面 ， 权 因子 起 到 i 
具有 对 称 性 ， 不 会 相差 太 大 ， 这 样 构造 出 来 的 曲面 不 会 很 扭曲 ， 应 该 是 对 称 的 . 
基于 这 种 情形 ， 下 面 考虑 一 种 特殊 的 有 到 


E Bézier 极 小 曲面 . 这 就 涉及 到 权 
周 整形 状 的 作用 ， 通 过 调和 求 出 的 权 因子 


E Bézier 极 小 曲面 设计 ， 那 么 应 怎样 求解 呢 ? 


是 否 可 以 利用 (3.2.4) 直 接 求 导 来 求解 呢 ? 下 一 节 将 对 此 加 以 阐述 . 


§3.3 有 限 差分 法 构造 有 理 Bézier 极 小 曲面 


直接 利用 Bézier 基 的 线性 无 关 性 来 寻找 控制 顶点 之 间 的 关系 ， 所 构造 的 曲面 是 较 精 而 


的 , 但 是 涉及 到 高 阶 非 线 性 方程 组 的 求解 ， 所 以 转换 思维 寻求 其 它 的 方法 来 解决 ， 采用 有 


限 差 分 法 . 虽然 上 只 是 一 种 逼近 而 不 是 精确 解 ,但 求解 过 程 是 求解 线性 方程 组 . 线性 方程 组 有 


其 特殊 优势 ， 求 解 快速 ， 在 实际 曲面 造型 


应 用 比较 多 . 


利用 有 限 差分 法 求解 这 一 问题 具体 思路 如 下 : 
将 所 求 的 调和 的 有 理 Bézier 曲面 划分 成 M x M 的 控制 网 格 ， 再 将 曲面 上 的 点 进行 离 


Mo WILK lly = G-D*1/M , v, e G-D*VM OBL j 21,2, M +L 


阶 偏 导 和 二 阶 偏 导 可 以 用 以 下 式 子 来 表示 : 


X ilj * i,j X hjd “ij 
, 3 
u h v h 
Y W Wij Y ija Wij 
, * 
d h ^ h 
X Xs 72% tX aug X Xin 2K tX; 
uu 2 > vv p , 
h h 
Y Wi 2w ij FW ji j W ijn 2w Lj TW 5j 
uu h? a yv h? 


由 偏 微 分 知识 可 知 : 24M ol, h—0 
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X; M p ELLO d 
+l, z JH : 
t 3 tJ > X ， iJ : LJ > X,, 
itl,j i,j i,j+l i,j 
h > 了， h >Y,, 
X ig Xi tX; X X pja TA tX, X 
h? > uu? h? > vv? 
Wij 2W /十 Y W i ja 2w; +W; ja Y 
h? > uu? h > vv’ 


(3.3.1) 


vo 


i=0 j=0 


tw, &u-(-D*1/M, v=(j-1)*1/M 时 ， 了 的 函数 值 


Iw, = 了 Gov) -Y Y B (G-0/M)B(G-D/M)w;. 


FMW, x; dtu =(G@-1*1/M, v-(o-D*I/M 时 ， 的 函数 值 . 


Mx, =X(u,v,) = 9. BG - 0 MOBIG - DI MDP,. 


i-0 j-0 

三 4H 

HR, + 及 ,=0 可 得 
w iau 0 i+l,j -2x ij +x ia) *W La (X i,j4l -2x ij +x Wa) 
—X ia, jW iW BS 2w QT MW aj) TH WG ja 
(Wij TZW FW xi 7X JW 44 — Wi) (3.3.2) 
Tu; TX POW ig 7Wi))Wwi;t2(wi;-w;;) 

' CURA 
HW 7W,) )*,; 70 
根据 第 二 章 的 知识 : 如 果 分 母 是 调和 的 ， 给 定 其 一 对 边界 的 权 因子 就 可 以 求 出 所 有 的 


RAF, ZA 


点 和 权 


Aly, ANH 


(3.3.2) 组 成 线性 方程 组 可 以 计算 4H 


解 是 唯一 的 . 最 后 将 所 求 出 的 点 连接 就 得 到 所 求 的 有 理 调和 Bézier 


fF 对 于 任意 给 定 的 U,v w ;都 是 已 知 的 , 给 定 有 下 


E Bézier 曲面 边界 的 控制 顶 


所 有 的 曲面 上 的 点 值 x，/w, ,并 且 此 


U y 


多 ， 其 曲面 越 光 请 ， 也 就 越 接近 于 有 
这 一 方法 ， 并 对 其 绝对 平均 曲率 图 进行 分 析 


Md 


里 Bézier 极 小 
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面 . 当 划 分 的 份 数 越 


将 通过 实例 来 验证 


H, 


面 . E p 
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第 四 章 有 理 Bkzier 极 小 曲面 求解 实例 


面 求解 实例 ， 并 对 其 绝对 平均 曲率 图 进行 综合 分 析 比 
， 其 平均 曲率 的 最 大 值 和 绝对 平均 曲率 的 平 


本 章 给 出 了 低 阶 有 理 Bézier 极 小 
较 ， 从 中 发 现 用 有 限 差 分 法 所 求 得 的 数值 解 
此 它们 都 可 以 看 作 是 近似 的 极 小 曲面 . 在 这 些 实例 中 ， 可 以 发 现 分 
面 . 


均值 都 比较 接近 于 零 ， 因 
母 调 和 的 情形 绝对 平均 曲率 更 接近 于 零 ， 更 加 逼近 于 有 到 


nj 
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$4.1 低 阶 求解 实例 


实例 1. 给 定 所 求 曲 面 的 控制 顶点 矩阵 (已 )， 如 下 : 


(Lb 25 0), t=]; 0, 2.5), (-L =2, 0) 
(0,2,2), | unknown,  (0,—2,2) 
(1, 2, 0), (1,0,2.5), (1, 2,0) 


ZDEBUS TREE (Wi) 5,3 如 下 : 


25 1, 2 
unknown, unknown, unknown |. 
2; 1, 2 
分 母 调和 的 情形 ; 
基于 第 二 章 中 的 算法 ， 利 用 分 母 调 和 所 求 得 的 权 因子 矩阵 如 下 : 
2: 1. 2 
3,°32,. 2 
2, lL 2 
当 曲 面 分 成 5x5 控制 网 格 , HI M = SIN, 需要 求解 16 个 线性 方程 , 得 到 曲面 上 的 25 个 点 ， 


(4.1.1). 


把 这 些 点 连接 起 来 就 得 到 离散 的 调和 的 有 理 Bézier 曲面 图 
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图 4.1.1 


当 曲 面 分 成 10x10 控 制 网 格 ， 即 WM =10 时 ， 需 要 求解 81 个 线性 方程 ， 得 到 曲面 上 的 121 


个 点 ， 把 这 些 点 连接 起 来 就 得 到 划分 10 份 的 离散 的 调和 的 有 理 Bézier 曲面 图 (4.1.2). 


4.1.2 


当 划 分 的 份 数 越 多 ， 所 求 的 曲面 越 光 清 . 为 了 从 客观 上 分 析 所 求 出 的 曲面 ， 就 要 分 析 
该 曲面 的 绝对 平均 曲率 图 ， 绝 对 平均 曲率 越 小 ， 所 求 出 的 曲面 也 就 越 接近 于 要 求 的 极 小 曲 


面 . 
1 LG-2MF + NE 
由 微分 几何 的 基本 知识 可 知 : 曲面 的 平均 aH SE Lm 


其 中 : EQFQG 表示 曲面 R 的 第 一 类 基本 量 ;，L,M,NN 表示 曲面 RR 的 第 二 类 基本 量 . 


> > > 


L=R,,:N.M=R,,: 


uu 


示 曲 面 R 的 法 线 量 ， 具 体 表达 式 如 下 : 


zZ 
Ku 
> 
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(RxR,) (R,xR,) 
RxR| VEG-F’ 
作 离 散曲 面 时 ， 是 通过 离散 曲面 求 出 曲面 上 点 的 值 ， 但 不 能 求 出 曲面 的 控制 顶点 ， 这 


N= 


也 是 造型 方法 的 缺陷 . 但 是 可 以 通过 拟 合 来 观察 其 曲率 的 总 体 变 化 趋势 ， 具 体 做 法 如 下 : 


对 于 二 次 有 理 Bézier 曲面 ， 未 知 控制 顶点 是 PP. 


取 所 求 出 的 曲面 上 的 一 个 非 边 界 上 的 点 值得 到 一 个 方程 来 求 出 了 ,. 如 曲面 分 成 


10x10 控 制 网 格 时 , 可 取 w = 0.5;v =0.5, 通过 上 述 算法 离散 曲面 ,R(0.5,0.5) 是 已 知 的 . 


2 2 B; 00)B) (vw, P; 


BE Ru, v) = EE — — À Wt M ea: 
B; (u) B; (v)w; 
i-0 j-0 
i=2 j-2 i-2 j=2 
B; (0.5) B (0.5)w,P, = R(0.5,0.5) B; (0.5) B7 (0.5)w, 
i-0 j=0 i-0 j=0 
经 整理 可 得 : 
i-2 j-2 2 2 
R(0.5,0.5) > 2 B; (0.5) B (0.5)w, — > >, B; (0.5) B5 (0.5)w, P, 
i-0 j=0 i-0 a 
2 
- 3, B? (0.5) BL (0.5)w, P, 
i=0 
P = idi A 4.1.1 
i B? (0.5)B? (0.5), ma: 


通过 (4.1.1) 求 出 近似 的 P| =[-0.0444 0.0584 -0.8756] .依据 平均 曲率 公式 ， 得 到 该 构 


造 曲面 的 平均 曲率 图 ， 从 中 观察 出 平均 曲率 的 绝对 变化 程度 如 图 4.1.3. 


4.1.3 
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其 它 情形 ( 权 因 子 任意 取 ): 
权 因 子 起 到 调整 曲面 形状 的 作用 ， 曲 面 会 向 权 因 子 较 大 的 控制 项 点 鼓 起 . 虽然 权 因 


之 间 相 差 不 是 很 大 ， 并 且 是 对 称 的 . 


反例 1.1 : 

所 取 的 权 因 子 和 矩阵 如 下 : 
2. 1, 2 
6, Ll 6| 
2, 1, 2 


当 曲 面 分 成 10x10 控 制 网 格 ， 即 M —10 时 ， 得 到 离散 的 曲面 图 4.1.4. 


图 4.1.4 


类 似 于 调和 情形 ， 求 出 控制 顶点 PB 的 近似 值 [-0.4090 0.8340 -4.6080] ， 再 求 得 


其 绝对 平均 曲率 图 4.1.5: 


图 4.1.5 
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反例 1.2 : 

所 取 的 权 因 子 和 矩阵 如 下 : 
2 
5, 3, 5. 
2, L 2 


当 曲 面 分 成 10x10 控 制 网 格 ， 即 WM = 10 时 ， 得 到 离散 的 曲面 图 4.1.6 . 


图 4.1.6 


类 似 于 调和 情形 ， 求 出 控制 顶点 近似 值 P, =[-0.1874 0.0653 -0.6445]， 


再 得 其 绝对 平均 曲率 图 4.1.7. 


图 4.1.7 


例 2:〈 三 阶 情形 ) 


给 定 所 求 曲面 的 控制 顶点 矩阵 (已 ) ,如 下 ， 
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(-1, 2, 0)， (-1, 0, 2.72), (-1,0, -2.72), (-1, -2, 0) 
(-0.5, 2, 2), | unknown, unknown, (-0.5, -2, 2) 
(0.5,2,2), | unknown, unknown, (0.5, -2, 0) | 

d,2,0),  (10,2.72, (10,2.72, (d, -2, 0) 


ERAT AREE QV), 如 下 : 


2, 1.5, 1.5, 2 
unknown, unknown, unknown, unknown 
unknown, unknown, unknown, unknown 


ps 1.5, 1.5, 2 
分 母 调和 情形 : 
利用 分 母 调和 求 得 的 权 因 子 和 矩阵 如 下 : 


2. 1.5, 1.5, 2 
25s 2. 2, 2.5 
2.5, 2, X Sr 
2, 1.5, 1.5, 2 


当 曲 面 分 成 10x10 控 制 网 格 , 即 M —10 时 ， 需 要 求解 81 个 线性 方程 ， 得 到 曲面 上 的 
121 个 点 ， 把 这 些 点 连接 起 来 就 得 到 划分 10 份 的 离散 的 调和 的 有 理 Bézier 曲面 图 4.1.8. 


图 4.1.8 
类 似 于 二 阶 情 形 ， 通 过 拟 合 来 观察 其 曲率 的 总 体 变化 . 取 所 求 的 曲面 上 的 非 边界 上 的 


四 个 点 的 值 组 成 方程 组 来 求 出 近似 的 P,P; Py, P. 


当 将 曲面 划分 为 3 份 时 ， 其 体 求解 过 程 如 下 : 
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Y Y B WB; ()w;P, 
E 得 一 个 线性 方程 组 : 
B? (u) Bi )w; 


t 


i=0 j=0 


i=3 j=3 


i , 


lud 1 
B B OWP; = RG 


i=0 j=0 


3 


1 2 
i 0 


i=0 j= 


i=3 j=3 


2 1 2 
B; (JB! (2w,P, = RE, 
5 CE 3 


t 


3 


ey) d dca cL 
Bi (B (Dw. P. = R(E, B? (2) B? (Sw, 
, iG) ; OW y G 3) : iG) iG) y 


i=3 j 
j=0 


i=0 j= i=0 j 


H ERIAREN ER lor DUI P, P, Pa Pp 值 如 下 : 


P,,=[-0.4904 0.3663 0.4298], 
P,,=[-0.4776 -0.2238 0.3813], 
P,,=[0.4580 0.3170 -0.0539], 
P,,=[0.4303 -0.2103 0.4835]. 


从 而 得 到 其 绝对 平均 曲率 图 4.1.9: 


图 4.1.9 
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其 他 情形 : 

反例 2.1 

所 取 的 权 因 子 和 矩阵 如 下 : 
2, 1.5, 1.5, 2 
5, 4 4 5 
5, 4, 4 5 
2, 1.5, 1.5, 2 


当 曲 面 分 成 $x$ 控制 网 格 ， 即 M =5 时， 需要 求解 16 个 线性 方程 ， 得 到 曲面 上 的 25 
个 点 ， 把 这 些 点 连接 起 来 就 得 到 划分 5 份 的 离散 的 调和 的 有 理 Bézier 曲面 图 4.1.10. 


图 4.1.10 


当 曲 面 分 成 10x10 控 制 网 格 ， 即 M =10 时 ， 需 要 求解 81 个 线性 方程 ， 得 到 曲面 上 的 
121 个 点 ， 把 这 些 点 连接 起 来 就 得 到 划分 10 份 的 离散 的 调和 的 有 理 Bézier 曲面 图 4.1.11. 


图 4.1.11 


求 得 其 未 知 控制 顶点 了 1,P, Pa Py IMEA F: 
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P=[-0.5631 0.3813 1.1635], 
P,=[ -0.6003 -0.2207 1.5255], 
P,,=[ 0.4405 0.2883 0.3882], 
P,,=[ 0.4085 -0.2004 0.2301]. 


进而 得 到 其 绝对 平均 曲率 图 4.1.12: 


图 4.1.12 


反例 2.2 
所 取 的 权 因 子 和 矩阵 如 下 : 


2. 1.5, 1.5, 2 
13. L 1, 1.5 
1.5, L 1 15[ 
2. 1.5, 1.5, 2 


当 曲 面 分 成 10x10 控 制 网 格 ， 即 M —10 时 ， 得 到 其 曲面 图 4.1.13 . 


图 4.1.13 
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求 得 其 未 知 控制 顶点 忆 1,P, P, Py 近似 值 如 下 : 


P,-[0.0736 -0.0297 -0.0172], 
P,-[0.2914 0.2901 -1.2228], 
P,-[-0.1180 0.0954 -1.1729], 
P,,=[ -0.4537 0.3990 2.4958]. 


进而 可 以 得 到 其 绝对 平均 曲率 图 4.1.14: 


实例 3 (四 阶 情形 ) 
给 定 所 求 曲面 的 控制 顶点 矩阵 ( 忆 Jog 如 下 : 


(0,0,0), (1,0,0), (4,00, (3,0,0), (40,0 | 
(0,1,0), unknown, unknown, unknown, (4, 1, 0.5) 
(0,2, 0), unknown, unknown, unknown, (4,2,1) 

(0,3,0), unknown, unknown, unknown, (4, 3, 0.5) 
(004,0, (14,05, (2,4,1), (,4,0.5), (4, 4, 0.2) | 


给 定 的 权 因 子 矩 阵 (QV), 如 下 : 


1.5, l, 1, 1, 1.5 
unknown, unknown, unknown, unknown, unknown 
unknown, unknown, unknown, unknown, unknown |. 
unknown, unknown, unknown, unknown, unknown 


1.5, 1, l, l, 1.5 ] 


调和 情形 : 
利用 分 母 调和 所 求 的 权 因子 矩阵 如 下 : 
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1.5, 1, l, l, 1.5 
2.5, 1.625 1.5 1.625 2.5 
3, 2, 1.8333, 2, 3 
2.5, 1.625, 1.5, 1.625, 2.5 
1.5, l, l, l, 1.5 


当 曲 面 分 成 10x10 探 制 网 格 ， 即 M =10 时 ， 需 要 求解 81 个 线性 方程 ， 得 到 曲面 上 的 
121 个 点 ， 把 这 些 点 连接 起 来 就 得 到 划分 10 份 的 离散 的 调和 的 有 理 Bézier 曲面 图 4.1.15. 


图 4.1.15 
类 似 于 三 阶 情形 ， 通 过 拟 合 来 观察 其 曲率 的 总 体 变化 . 取 所 求 出 的 曲面 上 的 非 边 界 上 


的 九 个 点 的 值 组 成 方程 组 来 求 出 近似 的 Pj, P. Pry, P,P,, P,P,P,P;. 


如 将 曲面 分 成 10 份 时 : 
u = 0.2,v = 0.3, R(0.2,0.3);u = 0.2, v = 0.5, R(0.2, 0.5); 
u =0.2,v = 0.7, R(0.2,0.7);u = 0.5, v = 0.2, R(0.5, 0.2); 
Eiu = 0.5, v = 0.2, R(0.5,0.2);u = 0.5, v = 0.5, R(0.5, 0.5); 
u = 0.5, v = 0.8, R(0.5,0.8);u = 0.8, v = 0.3, R(0.8, 0.3); 
u = 0.8, v = 0.7, R(0.8,0.7). 
i=4 j=4 
> >》B WB OwP, 
B Rv) == 可 得 九 个 线性 方程 : 


EN 
T 


i=4 j=4 


B*(0.2)B+(0.3)w,P, = R(0.2,0.3), 》 B (0.2) B: (0.3)w, 
0 =0 


i=4 j 
i j^ 


i-0 j i-0 j 
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求解 上 述 九 个 线性 方程 求 出 Pj, P, Pa, Pa Pros Pass 


i=4 j=4 


»' B^(0.5)B* (0.2), P, = R(0.5,0.2) 
=0 


j-4 j 


i-0 j 


进而 得 到 


4 


l 


i=4 j 
j=0 


i=0 j 
i=4 j=4 


t 
i-0 j-0 


t 


i=4 j-4 


了 


i=0 j=0 


iA j-4 
1 


i=0 j=0 


l 


( 


P,=[1.0000 
P,=[2.0000 
P,=[3.0000 
P,,=[1.0000 
P,,=[2.0000 
P,,=[3.0000 
P,,=[1.0000 
P,,=[ 2.0000 
P,,=[3.0000 


绝对 平均 曲率 图 4.1.16. 


1.0000 
1.0000 
1.0000 
2.0000 
2.0000 
2.0000 
3.0000 


i=4 j=4 


1 B; (0.2); (0.5)w;P; = R(0.2,0.5) 9 2 Bi (0.2) B; (0.5)w,. 


i-0 j=0 


i=4 j-4 


B; (0.2)85(0.7)w, P; = R(0.2,0.7) B; (0.2)B;(0.7)w,. 


l 


i=0 j=0 


=4 j=4 


B‘ (0.5) B (0.2)w;. 
0 


t 


ll 
`. 
I 


4 j=4 


B‘ (0.5) B (0.5)w, P, = 及 (0.5,0.5)》 》 B^ (0.5)B4(0.5)w;. 


t 
i-0 j=0 


i=4 j=4 


B; (0.5) B; (0.8)w,P, = R(0.5,0.8) 9 >, B‘ (0.5)B! (0.8)w;,. 


i-0 j=0 


i=4 j-4 i=4 j=4 
B; (0.8) B (0.3)w, P, = R(0.8,0.3) >) >) B‘ (0.8) B; (0.3)w,. 

-0 j-0 i-0 j-0 

i=4 j-4 i=4 j-4 
B; (0.8)B5(0.5)w, P, = R(0.8,0.5) B; (0.8)B(0.5)w,. 

=0 j-0 i-0 j-0 

4 i=4 j=4 


B; (0.8) Bi (0.7)w,P, = R(0.8,0.7) 9 > Bi (0.8) B; (0.7)w;. 


i-0 j-0 
P, P, Py, 如 下 : 


-0.4318], 

1.2335], 
-0.9739], 
0.5654], 
-0.5270], 
1.2616], 
0.0282], 


3.0000 0.5935], 


3.0000 
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4.1.16 
其 它 情形 : 
反例 3.1: 
所 取 的 权 因 子 和 矩阵 如 下 : 
15, l 1, 1, 15 
3, 2, 15, 2, 3 
4, 3, 2 3, 4 
3, 2, 15, 2, 3 
15, l 1, 1, 1.5 


当 曲 面 分 成 10x10 控 制 网 格 ， 即 M =10 时 ， 需 要 求解 81 个 线性 方程 ， 从 而 得 到 曲 
面 上 的 121 个 曲面 上 的 点 的 值 ， 可 得 到 离散 的 有 理 调和 Bézier 曲面 图 4.1.17 . 


图 4.1.17 


34 


浙江 大 学 硕士 学 位 论文 


P,-[1.0001 -0.0782 -0.4336], 
P,,=[0.9999 2.3504 0.5163], 
P,,=[1.0001 3.1551 0.1376], 
P,=[1.9998 0.3776 1.0735], 
P,,=[2.0003 1.9467 -0.5867], 
P,,=[1.9998 3.0074 0.8264], 
P,,=[3.0003 -0.0782 -1.0682], 
P,,=[2.9996 2.3504 1.4330], 
P,,=[3.0003 3.1551 0.0494]. 


从 而 得 到 其 曲面 的 绝对 平均 曲率 图 4.1.18: 


图 4.1.18 


这 中 只 给 出 了 所 求 出 的 曲面 图 ， 虽 然 给 出 了 其 对 应 的 绝对 平均 曲率 图 ， 但 没有 从 客 
观 上 加 以 分 析 . 


$ 4.2 曲率 图 的 综合 分 析 


首先 要 说 明 的 是 : 作 平 均 曲 率 图 时 , 将 曲面 分 成 21x21 的 网 格 , 得 到 曲面 上 的 441 
个 绝对 平均 曲率 的 值 . 

TE: 多 表示 未 知 的 权 因 子 ，M 表示 绝对 平均 曲率 的 最 大 值 ，A 表示 绝对 平均 曲率 的 
平均 值 ;已 表示 概率 ,可 以 定义 为 满足 条 件 的 曲率 面 上 的 点 的 个 数 比 上 总 的 个 数 441; 
F 表示 曲面 情况 ，7 表示 分 母 调 和 情形 ， 鳌 和 1 克 分 别 表示 其 他 两 种 权 因 子 任意 取 的 
Æ, H 表示 绝对 平均 曲率 . 
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第 一 组 实例 (二 阶 ): 

调和 情形 (7 P H 的 最 大 值 小 于 全 和 7 的 最 大 值 , 3f ELI P H PERS ELIT HI 
IULII PH. 
从 统计 的 角度 来 说 : 

对 于 情形 7， H < 0.05 WMH 22.68%, H < 0.1 的 概率 为 60.32%,, H < 0.5 
的 概率 为 85.49%， 万 > ] 的 概率 为 1.59% ; 


Merwe, H <0.05 的 概率 为 1.81%， 石 < 0.1 的 概率 为 10.20%, H <0.5 
的 概率 为 77.78%， 囊 >1 的 概率 为 3.40%; 


WEBI, H «0.05 WAN 11.11%, H < 0.1 的 概率 为 32.20%, H «0.5 
WIESE 76.87%, H » MESE 9.5296. 


A Ext rm wp EUR USE T SEECTEUE IE RU III SEU, OH «0.5 的 概率 就 有 


85.49%, F 22.68%) H < 0.05 EEFE, Ie TERI TEE 更 接近 于 要 求 的 有 理 Bézier 
调和 曲面 . 


表 ge 

F W M A P(H <0.05) | P(H <0.1) 
I [3, 2, 3] 1.1279 0.2299 100 266 
441 441 
II [6, 1, 6] 1.3458 0.3464 8 45 
441 441 
II] [5, 3, 5] 1.8460 0.3463 49 142 
441 441 
F P(H «0.5 P(H >=0.5) | P(H >=1) Ideal(H) 

7 377 Ei ES 
441 441 441 0 

1 343 98 15 
441 441 441 0 

m 339 102 a 
441 441 441 0 
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第 二 组 实例 (三 阶 ): 
调和 情 
SEF IL 


Zee) 


大 值 是 2.2319. 
从 统计 的 角度 来 说 : 
对 于 情形 了 ， 
(MEA 68.71%, H > 1 的 概率 为 3.85% ; 

而 对 于 
的 概率 为 60.09%， 万 > 1 的 概率 为 14.74%; 
对 于 情形 IM , 
(MEA 71.66%, H >] 的 概率 为 1.59% 
从 上 述 数 据 ! 


H «0.05 的 概率 为 9.75%， 


H 的 最 大 值 小 于 开 和 1 的 最 大 值 ， 并 且 相差 插 大 . IP H 的 平 
AS, DPH eA A 的 平均 值 相差 很 小 ， 但 是 鼠 中 豆 的 最 


H «0.05 的 概率 为 11.11%， H <0.1 的 概率 为 21.32%,, H < 0.5 


ÆI, H < 0.05 的 概率 为 8.62%， 态 < 0.1 的 概率 为 19.95%, H «0.5 


H < 0.1 的 概率 为 19.27%, H < 0.5 


可 以 看 出 情形 7 要 比 情形 代 明 显 的 好 ， 和 7 相 比 相差 不 大 ， 但 是 


IIl $$ H BERKIT, IPH «0.18IH «0.05 的 概率 都 比 1 下 大 ， 总 体 来 说 了 还 是 
要 好 一 些 的 . IPH < 0.5 的 概率 就 有 60.09%， 有 11.11% H < 0.05 WEF, 
We LI A TT 更 接近 于 要 求 的 有 理 Bézier 调和 曲面 . 
在 表 二 中 ， 玉 ,到 , 丙 可 以 表示 成 如 下 形式 : 
2.5, 2; 2; 23 3; 4, 4, 5 Ts dió 15 
W= ; W, = M 47 . 
2.552; 2; 25 5, 4, 4, 5 1, 1.5, 1, 1.5 
表 二 
F W M A | P(H <0.05) | P(H «0.) 
I W, 1.3894 0.3811 49 94 
441 441 
II W, 2.1195 0.5041 38 88 
441 441 
III W, 2.2319 0.3874 43 85 
441 441 
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F | P(H«05) | P(H>=0.5) | P(H >=) Ideal(H) 
I| — 303 138 am 
441 441 441 0 
n 265 176 65 
441 441 441 0 
J| 316 125 T 
441 441 441 0 
第 三 组 实例 (四 阶 ): 
调和 情形 (7 A H 的 最 大 值 小 于 二 中 玖 的 最 大 值 ， 并 且 相 差 挺 大. 1PH 的 平 
均值 小 于 本 中 的 ， 并 且 本 中 石 的 平均 值 基本 上 是 了 中 的 2 倍 ， 从 这 点 来 看 情形 工 是 
很 表 近 于 极 小 曲面 的 . 
从 统计 的 角度 来 说 : 


对 于 情形 7 ， 
的 概率 为 100%， 万 > ] 的 概率 为 0% ; 
而 对 于 情形 I , 
H «0.5 的 概率 为 95.24%， 万 > 1 的 概率 为 3.17; 
从 上 述 数据 


H «0.05 的 概率 为 69.84%, 


可 以 看 出 情形 了 要 比 情形 元 明 显 的 好 ， 


I} H «0.05 的 概率 就 有 


H «0.05 的 概率 为 78.68%， < 0.1 NEN 90.70%, H «0.5 


H «0.1 的 概率 为 82.09%, 


78.68%, H < 0.5 BE 7g 100%, H » LEGES 2g 0%, Mati n] AE I SÉ I SE 


接近 于 要 求 的 有 理 Bézier 调和 曲面 . 


在 表 三 中 ， 玉 ,到 , 丙 可 以 表示 成 如 下 形式 : 


2.5, 1.625, 1.5, 1.625, 2.5 3$. 2. dos 2; 

W,=| 3, 2, 1.8333, 2, 3,W-|43 2 3, 

2.5, 1.625, 1.5, 1.625, 2.5 3. 2. L5 2, 

表 三 : 
F W M A P(H «0.05) | P(H <0.1) 

I W, 0.4486 0.0382 347 400 
441 441 
II W, 1.1346 0.3123 308 362 
441 441 
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F P(H <0.5) | P(H>=0.5) | P(HS=1) Ideal(H) 
I 441 0 0 
441 441 441 0 
H 420 21 14 
441 441 441 0 
基于 上 述 三 个 表格 的 分 析 ， 可 知 分 母 调 和 情形 所 构造 出 来 有 理 Bézier 调和 曲面 其 绝对 


平均 曲率 变化 小 ， 更 通 近 于 极 小 曲面 . 在 上 述 的 造型 过 程 中 ， 为 了 保证 所 构造 曲面 的 凸 包 


性 质 并 避免 w(u,v) =0， 常 规定 wj > 0. 对 于 有 理 Bézier 调和 曲面 ， 权 因子 设计 有 着 很 重 


要 的 作用 .实验 表明 : w, 增 大 时 曲面 会 向 控制 顶点 靠 近 ， 从 而 调整 曲面 形状 


在 实际 的 造型 中 ， 为 了 保证 所 构造 的 曲面 比较 平坦 ， 要 求 权 因 子 要 匀称 ， 不 能 相差 太 


大 . 在 理论 上 ， 对 于 给 定 四 条 有 理 Bézier 边界 的 控制 顶点 及 第 一 排 和 最 后 一 排 权 因子 ， 利 


i=n j=n 
用 7 =X 》 B"(w)B" ()w, 调和 所 求 出 的 权 因子 分 布 比较 均匀 ， 相 差 不 会 很 大 ， 这 样 就 


i=0 j=0 


可 以 保证 所 求 出 的 曲面 不 会 很 扭曲 . 
对 于 上 述 反 例 ， 所 取 的 对 称 权 因 子 相差 不 是 很 大 ， 那 么 如 果 中 间 权 因子 取 的 是 对 称 但 


相差 很 大 时 造型 曲面 会 怎样 呢 ? 


4 


给 定 所 求 曲面 的 控制 顶点 矩阵 { 忆 ,} Au: 


| (1,2,0), (1 0, 2.72)， (-1 0, -2.72), (-1, -2,0) 
(-0.5, 2,2), | unknown, unknown, (-0.5, -2, 2) 
(0.5, 2, 2), unknown, unknown, (0.5, -2, 0) 
| (2,0)， (0,2.72)， (1,0, -2.72), (1, -2,0) 


4 


给 定 的 权 因子 矩阵 {W, 和 如下: 


2, 1.5, 1.5, 2 
unknown, unknown, unknown, unknown 
unknown, unknown, unknown, unknown | 

2; 1.5, 1.5, 2 


取 权 因子 矩阵 如 下 : 
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2: 2 
9 L 1 9 
9, 9| 
2. 1.5, 1.5, 2 

当 曲 面 分 成 10x10 控 制 网 格 , 即 M = 10 时 ， 得 到 离散 的 有 理 调和 Bézier 曲面 图 4.2.1 . 


图 4.2.1 


很 容易 看 出 其 曲面 不 是 很 平坦 ， 有 些 扭曲 ， 由 于 Wio,Woo,Wis, Wo EEK, HERS 


HŽ Po Po Pra. 忆 , 鼓 起 ， 使 得 在 这 儿 个 点 的 绝对 平均 曲率 变化 大 . 这 种 情形 下 即使 曲面 


大 部 分 绝对 平均 曲率 很 小 ,但 是 其 绝对 平均 曲率 的 最 大 值 会 很 大 . 

给 定 有 理 Bézier 曲面 控制 边界 ， 求 解 基 于 该 边界 的 面积 最 小 的 曲面 ， 需 要 求解 曲面 的 
E-LJjf£. 而 对 于 求解 其 一 工 方 程 有 很 多 方法 . 本 文中 采用 分 而 治之 的 方法 ， 先 使 得 分 
母 调和 ， 即 


(了 + 了 )=0 ， 这 时 (3.2.2) 中 一 (了 - Y, ) XY =0. 从 而 使 得 所 构造 的 有 理 Bézier 


曲面 更 接近 于 有 理 Bézier 极 小 曲面 . 权 因 子 任意 时 ，(Y,, +¥,) 不 是 恒 等 于 零 ， 使 得 多 构 


造 出 的 曲面 绝对 平均 曲率 变化 大 ， 这 一 点 可 以 由 育 


AT RS) SEP JI LAE HA. 


= 
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第 五 章 总 结 和 展望 


本 
验证 了 


的 是 曲 
成 果 ， 
构造 出 


这 一 点 
XX — FA 


我 认为 今后 的 工作 可 


第 


KB, 有 理 


数 最 小 
线性 问 


AAA 


ES 


文采 用 有 


该 算法 的 可 行 性 . 当 


限 差分 方法 成 功 地 解决 了 有 到 
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—> 


面 上 的 点 , 而 不 是 其 曲面 的 控 人 


E Bézier 调和 曲面 的 造型 问题 ， 
面 划 分 的 网 格 数 越 多 , 其 解 就 越 接近 了 
数 很 多 时 ， 虽 然 精 度 高 了 ， 但 计算 量 会 增 大 . 此 方法 适 
的 造型 问题 ， 其 求解 过 程 很 


虽然 这 种 方法 比较 简单 ， 但 是 基于 所 求 的 


E Bézier 调和 


TL 


面 进行 差分 离散 ， 求 解 线性 


通过 数值 试验 
六 极 小 曲面 . 当 划 分 的 份 
用 于 任意 阶 的 有 到 
才 单 ， 与 低 阶 类 似 ， 都 是 求解 一 个 线性 方程 组 


面 


方程 组 所 得 到 


着 顶 点 .另外 一 方面 ， 基 于 Plateau- Bézier 问题 的 研究 


本 文 着 重 探 索 了 分 母 调和 的 情形 . 通过 三 组 实例 综合 分 析 ， 可 知 分 母 调和 情形 下 所 


的 有 理 Bézier 调和 曲面 最 接近 于 极 小 曲面 . 本文 对 有 理 


加 以 说 明 ， 没 有 给 出 严谨 的 证 明 . 我 尝试 过 对 此 加 以 记 


是 以 后 的 研究 目标 之 一 


以 考虑 以 下 几 个 方面 : 


: 由 于 有 
极 小 
化 问题 . 


里 极 小 


面 问题 最 后 要 转换 为 一 个 线 怕 
因此 ,如 何 寻找 各 种 不 同上 


题 ,， 力 是 一 个 需要 着 重 碳 
: 本 


HAE] 


构造 现 


外 形 曲面 


AUI] Jj IR]. 


到 有 理 B 样 条 ， 即 将 有 


代 建筑 中 的 顶棚 、 网 染 、 
实际 拼接 也 是 离散 的 ， 


该 关注 


网 格 


fil, 4 


形 表 示 都 需 


上 | 有理 Be&zier 双 调包 
限 差 分 方法 推广 到 有 到 

B=: 从 本 质 上 说 ,有 限 差 分 方法 是 一 种 离散 数据 建 模 方法 . 就 目 
而 言 ， 越 来 越 多 的 几何 造型 和 图 


需 
我 们 所 使 


HERI, 


J DU 上 


时 数据 都 是 离散 


面 的 研究 . 
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因此 , 研究 离散 化 的 极 小 曲面 非常 必要 . 


分 曲面 等 离散 极 小 


Bézier 极 小 


面 的 求解 算法 


E 明 ， 但 是 还 没有 满意 的 结果 ， 关 于 


面 问题 的 高 度 非 线 性 , 解决 方案 一 般 要 采用 非 线 性 的 方法 . 具体 
FE 约 束 条 件 下 的 有 限 维 变量 的 非 线性 
适当 的 线性 化 策略 将 该 类 非 线性 问题 转化 为 某 些 


目标 函 


文 用 有 限 差分 方法 解决 有 理 Bézier 极 小 曲面 造型 时 , 理论 对 象 来 自 于 分 片 连续 
间 . 因此 ,可 以 考虑 将 此 方法 直接 推广 至 


[的 曲面 设计 ， 也 可 以 将 
EB 样 条 极 小 曲面 的 设计 中 . 


前 的 计算 机 表示 方法 
要 借助 于 离散 化 曲面 模型 . 当 使 用 有 限 差 分 方法 


的 ， 而 这 些 建筑 的 
Lp, SU 
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本 人 在 探索 有 理 Bézier 调和 曲面 造型 的 过 程 中 ， 最 先 想到 的 是 直接 利用 求 导 的 方法 ， 


求 导 后 ,将 所 求 的 有 理 多 项 式 升 阶 成 同 次 数 的 有 理 Bézier 曲面 ， 再 利用 Bézier 基 的 线性 无 


关 性 得 出 相应 的 结论 . 在 推导 的 过 程 中 ， 发 现 直 接 做 nxn 阶 的 情形 是 很 困难 的 . 这 时 开始 


尝试 低 阶 情形 . 对 于 二 次 情形 ， 我 做 了 比较 细致 的 推导 . 
PO NN 
2:2, BrG)B; vw, P, 
二 阶 有 理 Bézier 曲面 的 一 般 表示 形式 为 ， R=, | 
B; (u) B; (v)w, 
i-0 j=0 
记 其 分 子 ， 分 母 为 
X - Y BuB’ Q)w,P, 
i-0 j=0 
7 = B; (u) B; (v)w, 
i-0 j-0 
1 0 0 
利用 Bézier Jt 5; EZ ITE HL RAO AR: (B2,B2,B2)= (,u,w?)|-2 2 0 
1 -2 1 
将 分 子 分 母 分 别 转化 为 容 基 形式 : 
i=2 j=2 : 
Xo u'v'a, 
i=0 j=0 
i=2 j=2 
Y = uv’b, 
i=0 j=0 
| am Ay, 0, 1 0 0| Wooo woPy Woo |Z -2 I 


其 中 | an a, 45,|2|-2 2 0| w,P, wiP: wp ||0 2 -2 
[Q5 05 Ay 1 -2 1| WP wb w5jP,]0 0 1 


| by, ba bo 1.0 0 | woo Wa Wo|1 -2 1 
0 Wie w wall0 2 -2 


|b, b b 1 -2 Tw wy wll0 0 1 


SW 
ST 
SW 
N 
Il 
LI 
N 
N 


(X,, +X,,¥? -XYY «Y, ) -2Y (X,Y,+ XY,)+2X(¥7+¥7)=0 (9 


将 分 子 表示 成 矩阵 形式 : 
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n 


对 于 分 母 : 
Do 
Y =|1 u u^ | b, 
b, 


X, 2 X, ? X,, * X, ? X, 2 Y, 2 Y, ^" uu? vv? uv 


Y, ,了 ,了 ,也 有 类 


似 矩 阵 表示 . 


利用 上 述 (*) 式 得 到 一 个 方程 ， 易 知 (*) 式 ! 
u,v 的 各 种 组 合 如 下 : 


u 的 最 高 次 数 为 6， v 的 最 高 次 数 也 为 6. cm 


L u uw uw uw w u‘ | [7 vy y? y? y* y) v‘ |= 
I v v v” v y v | 
u uv uv? uv? uv uv uv 
u^ uv wy? wy wyv* wy wy? 
u? wyv wy wy wy wy wy? (**) 
u^ uty uv? uv? uty? uv? wy? 
w wy wv wy wy wy wy? 
[us uv uv? uv? uv uv? uv 
有 理 Bezier HH MV ACE EEC AE 0 AF REI ZR PE DCSE TN, E) 
恒 等 于 0 UTE EA ALC) p HE IH MLE, 
从 而 得 到 方程 如 下 : 


h00 =-4*x10*(-2+2*w10)-4*x10-20 


h01 =—4* x10* (-4—4* wl0+4* wl 1) -8* x10 +32 —8* x11—8* x10*(—2+2* wl0) 
-(-8* x10—16-- 8* x11) *(-2+ 2*w10) +16* w104-8* (-2 4- 2* w10) 


h03 2 8*(-4—4* wl0+ 4* wl 1)? +16* x10--56—16* x12—(-8* x10—16 - 8* x11) 


*(442*wl0—4*wll+2* wl2)+16* wll+(— 


32 +32* wl0)*(4+2* wl0—4* wll+2* wl2) 


+8 * x10* (4-4* wl10+4* wll) +16* wl2—-(4* x10+8—8* x11+4* x12)* 


(4-4* wl0+4* wll) +32* x11—(8* x10+1 
—(—16* x10—324+16* x11) *(-4—4* wl0+4 


6—16* x11+8* x12) *(—2+2* wl0)—80* w10 
*wl1)-(16-16*w10)*(—4—4*wl10-- 4*wl1) 


42*(-8* x10—16-8* x11) *(-2-2*w10) -8* x10*(4-2*wl10—4*wl1-2*w12) 


类 似 地 ， 可 得 所 有 的 系数 阵 如 下 : 
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hy, h, hy h, i bs, hg; 


20> I1» I5 fhg» Mog» fhs» Myo 3 
hoo, bai, has, Pss gy» Mos, Mag’ 
yos hyss hs, a hy, hs, Mag: 
hso, Asi» Asz; hss, hsa, hss, hse; 


ho > h he, > has > ha > Nes d hs d 


1? 


c 
> 


1? 


二 S osse Maja hss, Mags ho, hos, Pee HE, h; RIRE E uv! 的 系数 . 


Iii RC, v) 调和 < hh==0(i,j — 0,1,......6). 


这 样 就 会 产生 40 个 非 线 性 方程 所 组 成 的 方程 组 . 别 说 是 求解 , 即使 判断 这 个 方程 组 是 
否 有 解 都 是 很 困难 的 . 
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